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Abstract. We consider minimal surfaces M which are complete, embedded 
and have finite total curvature in R 3 , and bounded, entire solutions with finite 
Morse index of the Allen-Cahn equation Am + f(u) = in R 3 . Here / = —W 
with W bistable and balanced, for instance W(u) = j(l — u 2 ) 2 . We assume 
that M has m > 2 ends, and additionally that M is non-degenerate, in the 
sense that its bounded Jacobi fields are all originated from rigid motions (this is 
known for instance for a Catenoid and for the Costa-Hoffman-Meeks surface of 
any genus). We prove that for any small a > 0, the Allen-Cahn equation has a 
family of bounded solutions depending on m — 1 parameters distinct from rigid 
motions, whose level sets are embedded surfaces lying close to the blown-up 
surface M a := « -1 M, with ends possibly diverging logarithmically from M a . 
We prove that these solutions are L°° -non- degenerate up to rigid motions, and 
find that their Morse index coincides with the index of the minimal surface. 
Our construction suggests parallels of De Giorgi conjecture for general bounded 
solutions of finite Morse index. 
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1. Introduction and main results 

1.1. The Allen-Cahn equation and minimal surfaces. The Allen-Cahn equa- 
tion in WL N is the semilinear elliptic problem 

Au + f{u) = inl", (1.1) 

where f(s) = —W'(s) and W is a "double-well potential", bi-stable and balanced, 
namely 

W(s) > if s ^ 1,-1, W(l) = = W(-l), W"(±l) = /'(±1) =: o\ > 0. 

" (1.2) 

A typical example of such a nonlinearity is 

f{u) = {l-u 2 )u for W(u) = -(1 -u 2 ) 2 , (1.3) 

while we will not make use of the special symmetries enjoyed by this example. 

Equation (|1.1|) is a prototype for the continuous modeling of phase transition 
phenomena. Let us consider the energy in a subregion region of 

Uv)= f ^\Vv\ 2 + ±.W(v), 

whose Euler-Lagrange equation is a scaled version of (jl.ip , 

a 2 Av + f(v) = in Q. (1.4) 

We observe that the constant functions u — ±1 minimize J a . They are idealized 
as two stable phases of a material in O. It is of interest to analyze stationary 
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configurations in which the two phases coexist. Given any subset A of fi, any 
discontinuous function of the form 

v* = xa - Xn\A (1.5) 

minimizes the second term in J e . The introduction of the gradient term in J a 
makes an a-regularization of u* a test function for which the energy gets bounded 
and proportional to the surface area of the interface M = dA, so that in addition to 
minimizing approximately the second term, stationary configurations should also 
select asymptotically interfaces M that are stationary for surface area, namely 
(generalized) minimal surfaces. This intuition on the Allen-Cahn equation gave 
important impulse to the calculus of variations, motivating the development of the 
theory of T-convergence in the 1970's. Modica [27] proved that a family of local 
minimizers u a of J a with uniformly bounded energy must converge in suitable sense 
to a function of the form (|1.5|) where <9A minimizes perimeter. Thus, intuitively, for 
each given A S (—1,1), the level sets [v a = A], collapse as a — > onto the interface 
dA. Similar result holds for critical points not necessarily minimizers, see [23]. For 
minimizers this convergence is known in very strong sense, see [2J[3J. 

If, on the other hand, we take such a critical point u a and scale it around an 
interior point G Q, setting u a (x) = v a (ax), then u a satisfies equation in an 
expanding domain, 

Aw a + f(u a ) = in a fi 
so that letting formally a — > we end up with equation (jl.lj) in entire space. The 
"interface" for u a should thus be around the (asymptotically flat) minimal surface 
M a = a _1 M. Modica's result is based on the intuition that if M happens to be 
a smooth surface, then the transition from the equilibria -1 to 1 of « Q along the 
normal direction should take place in the approximate form u a (x) ps w(z) where z 
designates the normal coordinate to M a . Then w should solve the ODE problem 

w" + f{w) =0 in K, w(-oo) = -1, w(+oo) = 1 . (1.6) 
This solution indeed exists thanks to assumption |T72]). It is strictly increasing and 
unique up to constant translations. We fix in what follows the unique w for which 

/ tw'(tfdt = . (1.7) 
Jr 

For example (|1.3[) . we have w(t) = tanh (i/v^)- In general w approaches its limits 
at exponential rates, 

w(t)-±l = 0(e- ff±|t| ) asi-*±oo. 

Observe then that 

J a (u a ) fa Area (M) [ \^-w' 2 + W(w)} 

which is what makes it plausible that M is critical for area, namely a minimal 
surface. 

The above considerations led E. De Giorgi [3] to formulate in 1978 a celebrated 
conjecture on the Allen-Cahn equation (jl.ip . parallel to Bernstein's theorem for 
minimal surfaces: The level sets [u — A] of a bounded entire solution u to (|1 . 1)1 . 
which is also monotone in one direction, must be hyperplanes, at least for dimension 
N < 8. Equivalently, up to a translation and a rotation, u = w(x\). This conjecture 
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has been proven in dimensions N = 2 by Ghoussoub and Gui [13], N — 3 by 
Ambrosio and Cabre pQ , and under a mild additional assumption by Savin |34j . A 
counterexample was recently built for N > 9 in [TOl E] , see also [21] • See [T2] 
for a recent survey on the state of the art of this question. 

The assumption of monotonicity in one direction for the solution u in De Giorgi 
conjecture implies a form of stability locally minimizing character for u when com- 
pactly supported perturbations are considered in the energy. Indeed, if Z = d XN u > 
0, then the linearized operator L = A + f'(u), satisfies maximum principle. This 
implies stability of u, in the sense that its associated quadratic form, namely the 
second variation of the corresponding energy, 

:= / |V</f - /»</> 2 (1.8) 

satisfies Q(i/j,ip) > for all ip ^ smooth and compactly supported. Stability is a 
basic ingredient in the proof of the conjecture dimensions 2,3 in [T| I13j. based on 
finding a control at infinity of the growth of the Dirichlet integral. In dimension 
N = 3 it turns out that 

/ |Vu| 2 = 0(R 2 ) (1.9) 

which intuitively means that the embedded level surfaces [u = A] must have a finite 
number of components outside a large ball, which are all "asymptotically flat" . 
The question whether stability alone suffices for property (II. 9|) remains open. More 
generally, it is believed that this property is equivalent to finite Morse index of 
the solution u (which means essentially that u is stable outside a bounded set). 
The Morse index m(u) is defined as the maximal dimension of a vector space E of 
compactly supported functions such that 

Q(V>,V)<0 for all ^e£\{0}. 

Rather surprisingly, basically no examples of finite Morse index entire solutions 
of the Allen-Cahn equation seem known in dimension N = 3. Great progress has 
been achieved in the last decades, both in the theory of semilinear elliptic PDE like 
(jl.ip and in minimal surface theory in K 3 . While this link traces back to the very 
origins of the study of (jl.ip as discussed above, it has only been partially explored 
in producing new solutions. 

In this paper we construct a new class of entire solutions to the Allen-Cahn 
equation in R 3 which have the characteristic (|1.9[) . and also finite Morse index, 
whose level sets resemble a large dilation of a given complete, embedded minimal 
surface M, asymptotically flat in the sense that it has finite total curvature, namely 

\K\dV < +oo 

where K denotes Gauss curvature of the manifold, which is also non- degenerate in 
a sense that we will make precise below. 

As pointed out by Dancer [7], Morse index is a natural element to attempt 
classification of solutions of (jl.ip . Beyond De Giorgi conjecture, classifying solutions 
with given Morse index should be a natural step towards the understanding of the 
structure of the bounded solutions of ([l.ip . Our main results show that, unlike the 
stable case, the structure of the set of solutions with finite Morse index is highly 
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complex. On the other hand, we believe that our construction contains germs 
of generality, providing elements to extrapolate what may be true in general, in 
analogy with classification of embedded minimal surfaces We elaborate on these 
issues in Cl2l 



1.2. Embedded minimal surfaces of finite total curvature. The theory of 
embedded, minimal surfaces of finite total curvature in R 3 , has reached a notable 
development in the last 25 years. For more than a century, only two examples of 
such surfaces were known: the plane and the catenoid. The first nontrivial example 
was found in 1981 by C. Costa, [U[S]. The Costa surface is a genus one minimal 
surface, complete and properly embedded, which outside a large ball has exactly 
three components (its ends), two of which are asymptotically catenoids with the 
same axis and opposite directions, the third one asymptotic to a plane perpendicular 
to that axis. The complete proof of embeddedness is due to Hoffman and Meeks 
[18] . In [19j [21] these authors generalized notably Costa's example by exhibiting 
a class of three-end, embedded minimal surface, with the same look as Costa's far 
away, but with an array of tunnels that provides arbitrary genus k > 1. This is 
known as the Costa-Hoffman-Meeks surface with genus k. 

Many other examples of multiple-end embedded minimal surfaces have been 
found since, see for instance [25] [36] and references therein. In general all these 
surfaces look like parallel planes, slightly perturbed at their ends by asymptotically 
logarithmic corrections with a certain number of catenoidal links connecting their 
adjacent sheets. In reality this intuitive picture is not a coincidence. Using the 
Eneper-Weierstrass representation, Osserman [21] established that any embedded, 
complete minimal surface with finite total curvature can be described by a conformal 
diffeomorphism of a compact surface (actually of a Riemann surface), with a finite 
number of its points removed. These points correspond to the ends. Moreover, after 
a convenient rotation, the ends are asymptotically all either catenoids or plane, all 
of them with parallel axes, see Schoen [35]. The topology of the surface is thus 
characterized by the genus of the compact surface and the number of ends, having 
therefore "finite topology" . 

1.3. Main results. In what follows M designates a complete, embedded minimal 
surface in R 3 with finite total curvature (to which below we will make a further 
nondegeneracy assumption). As pointed out in |22| . M is orientable and the set 
R 3 \ M has exactly two components S+, S-. In what follows we fix a continuous 
choice of unit normal field v(y), which conventionally we take it to point towards 
S+. 

For x = (xx,X2,Xz) — (x',xg,) € R 3 , we denote 



After a suitable rotation of the coordinate axes, outside the infinite cylinder r < Rq 
with sufficiently large radius Rq, then M decomposes into a finite number m of 
unbounded components Mi,...,M m , its ends. From a result in [35], we know 
that asymptotically each end of Mk either resembles a plane or a catenoid. More 
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precisely, M k can be represented as the graph of a function F k of the first two 
variables, 

M k = { y e R 3 / r(y) > 1%, Vs = F k (y') } 
where F k is a smooth function which can be expanded as 

F k (y') = a k log r + b k + b lk ^ + 0(r~ 3 ) as^+(», (1.10) 

for certain constants a k , b k , bi k , and this relation can also be differentiated. Here 

rn 

a\ < a 2 < . . . < a m , ^ a k = 0. (1-11) 

fc=i 

The direction of the normal vector v{y) for large r(y) approaches on the ends that 
of the X3 axis, with alternate signs. We use the convention that for r(y) large we 
have 

(-l) h 

v(y) = , =(V^fcfa / ),-l) ifyeM fe . (1.12) 

Let us consider the Jacobi operator of M 

J{h) A M h+\A\ 2 h (1.13) 

where \A\ 2 — —2K is the Euclidean norm of the second fundamental form of M . J 
is the linearization of the mean curvature operator with respect to perturbations of 
M measured along its normal direction. A smooth function z(y) defined on M is 
called a Jacobi field if 3{z) = 0. Rigid motions of the surface induce naturally some 
bounded Jacobi fields: Associated to respectively translations along coordinates 
axes and rotation around the X3-axis, are the functions 

z i(y) = Hv) ■ e i> yeM, i = 1,2,3, 

*4(tf) = (-l&, Ift.O) yeM. (1.14) 
We assume that M is non- degenerate in the sense that these functions are actu- 
ally all the bounded Jacobi fields, namely 

{ z€L°°(M) / J(z) = Q} = span{z 1 ,z 2 ,Z3,z 4 }. (1.15) 

We denote in what follows by J the dimension (< 4) of the above vector space. 

This assumption, expected to be generic for this class of surfaces, is known in 
some important cases, most notably the catenoid and the Costa-Hoffmann-Meeks 
surface which is an example of a three ended M whose genus may be of any order. 
See Nayatani [29, 30] and Morabito [28]. Note that for a catenoid, Z04 = so that 
J = 3. Non-degeneracy has been used as a tool to build new minimal surfaces for 
instance in Hauswirth and Pacard [17] , and in Perez and Ros [33] . It is also the basic 
element, in a compact-manifold version, to build solutions to the small-parameter 
Allen-Calm equation in Pacard and Ritore [52] . 

In this paper we will construct a solution to the Allen Cahn equation whose 
zero level sets look like a large dilation of the surface M, with ends perturbed 
logarithmically. Let us consider a large dilation of M, 

M a := oT x M. 

This dilated minimal surface has ends parameterized as 

M k>a = {yeR 3 I r(ay) > R , y 3 = a^F^ay') } . 
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Let (3 be a vector of given m real numbers with 

m 

= (/3 x ,...,/3 m ), X> = 0. (1.16) 

i=l 

Our first result asserts the existence of a solution u = u a defined for all sufficiently 
small a > such that given A € (—1,1), its level set [u a — A] defines an embedded 
surface lying at a uniformly bounded distance in a from the surface M a , for points 
with r(ay) — O(l), while its A:-th end, k = 1,. .. ,m, lies at a uniformly bounded 
distance from the graph 

r(ay) > R , y 3 = a' 1 F k (ay') + p k log \ay'\ . (1.17) 

The parameters (3 must satisfy an additional constraint. It is clear that if two 
ends are parallel, say a k +i = a k , we need at least that P k +\ — Pk > 0, for otherwise 
the ends would eventually intersect. Our further condition on these numbers is that 
these ends in fact diverge at a sufficiently fast rate. We require 

Pk+i - Pk > 4 maxlcrl 1 ,^ 1 } if a k+ % = a k ■ (1-18) 

Let us consider the smooth map 

X(y,z) = y + zv(ay), (y,t) £ M a xl (1.19) 

x = X(y,z) defines coordinates inside the image of any region where the map is 
one-to-one. In particular, let us consider a function p(y) with 

p(y) = {-l) k fi k log \ay'\ + 0(1), k = 1, . . . , m, 

and P satisfying Pk+i — Pk > 7 > for all k with a k = a k+ \. Then the map X is 
one-to-one for all small a in the region of points (y, z) with 

\z~q(y)\ < -+7log(l + |m/|) 
a 

provided that S > is chosen sufficiently small. 

Theorem 1. Let N = 3 and M be a minimal surface embedded, complete with 
finite total curvature which is nondegenerate. Then, given P satisfying relations 
(| 1 . 1 6\ and _?<$[) , there exists a bounded solution u a of equation (jj.jp . defined for 
all sufficiently small a, such that 

Ua(x) — w(z — q(y)) + Ota) for all x = y + zv(ay), \z — q(y)\<—, (1-20) 

a 

where the function q satisfies 

q(y) = {-l) k Pk log \ay'\ + 0(1) y € M k , a , k = 1, . . . , m. 

In particular, for each given A S (—1,1), the level set [u a — A] is an embedded 
surface that decomposes for all sufficiently small a into m disjoint components 
(ends) outside a bounded set. The k-th end lies at 0(1) distance from the graph 

V3 = a^ 1 F k {ay) + p k \og\ay'\. 

The solution predicted by this theorem depends, for fixed a, on m parameters. 
Taking into account the constraint Y^JLi Pi = this gives m — 1 independent 
parameters corresponding to logarithmic twisting of the ends of the level sets. Let 
us observe that consistently, the combination P £ Span{(ai, . . . , a m )} can be set 
in correspondence with moving a itself, namely with a dilation parameter of the 
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surface. We are thus left with m — 2 parameters for the solution in addition to 
a. Thus, besides the trivial rigid motions of the solution, translation along the 
coordinates axes, and rotation about the X3 axis, this family of solutions depends 
exactly on m — 1 "independent" parameters. Part of the conclusion of our second 
result is that the bounded kernel of the linearization of equation (jl.ip about one of 
these solutions is made up exactly of the generators of the rigid motions, so that in 
some sense the solutions found are L°°-isolated, and the set of bounded solutions 
nearby is actually m — 1 + J-dimensional. A result parallel to this one, in which 
the moduli space of the minimal surface M is described by a similar number of 
parameters, is found in [33] . 

Next we discuss the connection of the Morse index of the solutions of Theorem 
[T] and the index of the minimal surface M, i(M), which has a similar definition 
relative to the quadratic form for the Jacobi operator: The number i(M) is the 
largest dimension for a vector spaced E of compactly supported smooth functions 
in M with 

f \Vk\ 2 dV-[ \A\ 2 k 2 dV <0 for all keE\{0}. 
Jm Jm 

We point out that for complete, embedded surfaces, finite index is equivalent to 
finite total curvature, see [16] and also §7 of [22] and references therein. Thus, 
for our surface M, i(M) is indeed finite. Moreover, in the Costa-Hoffmann-Meeks 
surface it is known that i(M) = 21 — 1 where I is the genus of M. See [29], [30] and 

EE]. 

Our second result is that the Morse index and non-degeneracy of M are trans- 
mitted into the linearization of equation (II. ip . 

Theorem 2. Let u a the solution of problem (j_Z.jp given by Theorem^ Then for 
all sufficiently small a we have 

m(u a ) = i(M). 

Besides, the solution is non-degenerate, in the sense that any bounded solution of 

A0 + /' {u a )<j> = mR 3 
must be a linear combination of the functions Zi, i = 1, 2, 3, 4 defined as 
Zi = diU a , i = l,2, 3, Z 4 = -x 2 diu a + xid 2 u a . 

We will devote the rest of this paper to the proofs of Theorems 1 and 2. 
2. The Laplacian near M a 

2.1. The Laplace-Betrami Operator of M a . Let D be the set 

D = {y e M 2 I |y| > R }- 

We can parameterize the end of M as 

y € -D 1 — > y := Y k (y) = + F fc (y)e 3 . (2.1) 

and Fk is the function in (jl.lOp . In other words, for y = (y' , y 3 ) G Mk the coodinate 
y is just defined as y = y' . We want to represent AM~the Laplace-Beltrami operator 
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of M-with respect to these coordinates. For the coefficients of the metric gij on 
M k we have 

d Yl Y k = e I + 0(r- 1 )e 3 

so that 

fltf (y) = (9iYk, djY k ) = 6 tJ + O (r- 2 ) , (2.2) 
where r — |y|. The above relations "can be differentiated" in the sense that differ- 
entiation makes the terms 0(r~ J ) gain corresponding negative powers of r. Then 
we find the representation 



A M = t= d l (y/detg~g^d J )=A y + 0(r- 2 )d lJ +0(r- 3 )d l on M k . (2.3) 
Vdetgy 

The normal vector to M at y G Affc fc = 1, . . . , to, corresponds to 

= (-i) fc ^_L_= (^(yjei - e 3 ) , y = F fe (y) G M fe 
VI + |VF fe (y)| 2 

so that 

v{y) = {-I ) k e 3 + a k r- 2 + 0( r - 2 ) , y = Y k (y) G M k . (2.4) 

Let us observe for later reference that since div = 0(r~ 2 ), then the principal 
curvatures of M, k±, ki satisfy ki = 0{r~ 2 ). In particular, we have that 

\A(y)\ 2 = k 2 1 +k 2 2 =0(r- i ). (2.5) 

To describe the entire manifold M we consider a finite number N > m + 1 of local 
parametrizations 

yeWtC K 2 i — ► y = F fe (y), Y k G C°°(W fc ), = 1, . . . , N. (2.6) 

For fc = f, . . . , to we choose them to be those in (|2.ip . with = £), so that 
Y k (U k ) = M k , and % is bounded for fc = to + 1, . . . , N. We require then that 

N 

M=\J Y k {U k ). 

fc=i 

We remark that the Weierstrass representation of M implies that we can actually 
take N = m + 1 , namely only one extra parametrization is needed to describe the 
bounded complement of the ends in M. We will not use this fact. In general, we 
represent for y G Y k {U k ), 

A M = a ij (y)d ij + b° i (y)d i , y = Y k (y), yeU k , (2.7) 

Hi 

coefficients. For fc = 1, . . . , to we have 

a%{y) = 5 ij + 0{r- 2 ), b^O^ 3 ), as r(y) = |y| oo. (2. 



where is a uniformly elliptic matrix and the index fc is not made explicit in the 



The parametrizations set up above induce naturally a description of the expanded 
manifold M a = a~~ 1 M as follows. Let us consider the functions 

Y ka : U ka := oT x U k -> M a , y ^ Y ka {y) := a^Y^ay), k = 1, . . . , N. (2.9) 

Obviously we have 

N 

M a = |J Y ka {U ka ). 
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The computations above lead to the following representation for the operator A m q : 
Am„ = a° ij (ay)d ij + b° i (ay)d i , y = Y ka (y), yeU ka , (2.10) 
where a°- , 6? are the functions in (|2.7j) , so that for A; = 1, . . . , m we have 

a%=5 rj +0{r- 2 ), 6° = 0(r" 3 ), as r a (y) := |«y| -> oo. (2.11) 



2.2. The Euclidean Laplacian near M a . We will describe in coordinates relative 
to M a the Euclidean Laplacian A x , x £ R 3 , in a setting needed for the proof of our 
main results. Let us consider a smooth function h : M — > R, and the smooth map 
Xh defined as 

X h :M a xR^R 3 , (y,t)^X h (y,t):=y+(t + h(ay))v{ay) (2.12) 

where v is the unit normal vector to M. Let us consider an open subset O of M a x R 
and assume that the map Xh\o is one to one, and that it defines a diffeomorphism 
onto its image Af = Xh(0). Then 

x = X h (y,t), (y,t)eO, 

defines smooth coordinates to describe the open set Af in R 3 . Moreover, the maps 

x = X h (Y ka (y),t), (y,t) e {U ka x R) nO, k = l,...,N, 

define local coordinates (y, t) to describe the region Af. We shall assume in addition 
that for certain small number S > 0, we have 

O C {(y, t) I \t + h(ay)\ < - log(2 + r a (y) ) }. (2.13) 

a 

We have the validity of the following expression for the Euclidean Laplacian 
operator in Af. 

Lemma 2.1. For x — Xh(y,t), (y,t) S O with y = Y ka (y), y G U ka , we have the 
validity of the identity 

A x = d tt + A Ma - a 2 [{t + h)\A\ 2 + A M h]d t - 2aa% djhBu + 
a(t + h) [a^dtj - 2a ajj dihd jt + ab\ (d t - adihd t ) ] + 

a 3 (t + h) 2 b\d t + a 2 [a° tj + a(t + fe)ajj }d i hd j hd tt . (2.14) 

Here, in agreement with \2. 10$ , Am = a° i j(ay)dij + b° i (ay)di. 
The functions a\j , b\, b\ in the above expressions appear evaluated at the pair 
(ay , a(t + h(ay)) , while the functions h, dih, A^jh, \A\ 2 , a^, &° are evaluated at 
ay In addition, for k = 1, . . . , m, I = 0, 1, 

<4 = k 3 5 M + 0(r- 2 ), b\ = 0{r- 3 ), b\ = 0(r- 6 ) , 

as r a{y) — \ocy\ — * oo, uniformly in their second variables. The notation djh refers 
to dj[h o Y k ]. 
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We postpone proof of this fact for the appendix. The proof actually yields that 
the coefficients a}j and b\ can be further expanded as follows: 

a \j = a l 3 i a V, °) + «(* + h ) a if {ay, a{t + h)) =: a\f + a{t + h)a 2 3 , 
with <iy = 0(r~ 3 ), and similarly 

b) = b)(ay, 0) + a(t + h) bf\ay, a(t + h)) =: b)'° + a(t + h)b 2 , 

with bj = 0(r~ 4 ). As an example of the previous formula, let us compute the 
Laplacian of a function that separates variables t and y, that will be useful in Sj3] 
and §11. 

Lemma 2.2. Let v(x) — k{y) ip(t) . Then the following holds. 

A x v = kip" +ipA Ma k - a 2 [{t + h)\A\ 2 + A M h] kip' - 2aa° ij d j hd i kip' + 

a{t + h) [alfdijktp - 2a a]f djhdik ip' + a{b\ () d l k i/> - ab 1 l °d l h ki//)] + 

a 2 (t + h) 2 [dijdijkip - 2a o?- djhdik ip' + a{b\dik ip - abjd^ kip 1 )] + 

a^{t + h) 2 b\kiP' + a 2 [a° i3 + a(t + ^al^d.hdjhkiP" . (2.15) 



3. Approximation of the solution and preliminary discussion 

3.1. Approximation of order zero and its projection. Let us consider a func- 
tion h and sets O and M as in i j2.21 Let x = Xh(y, t) be the coordinates introduced 
in (|2.12[) . At this point we shall make a more precise assumption about the function 
h. We need the following preliminary result whose proof we postpone for £16.21 

We consider a fixed m-tuple of real numbers (3 = . . . , f3 m ) such that 

m 

Lemma 3.1. Given any real numbers /3i,...,/3 TO satisfying \3.1\ , there exists a 
smooth function ho(y) defined on M such that 

J (h ) = A M h + \A\ 2 h = inM, 

ho(y) = (—iyPj logr + 9 as r — > oo in Mj for all y € Mj , 
where 6 satisfies 

H0IU + Hr^lU < +oo . (3.2) 

We fix a function ho as in the above lemma and consider a function h in the 
form 

h = ho + hi . 
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We allow hi to be a parameter which we will adjust. For now we will assume that 
for a certain constant K, we have 

l|/ii|U-(M) + ||(l + r 2 )£>/ii|U=o (M) < Ka. (3.3) 

We want to find a solution to 

S(u) :=A x u + f(u) = 0. 

We consider in the region Af the approximation 

u (x) := w(t) = w{z - h (ay) ~ hi(ay)) 

where z designates the normal coordinate to M a . Thus, whenever (3j ^ 0, the level 
sets [uo = A] for a fixed A e (—1,1) departs logarithmically from the end a~ 1 Alj 
being still asymptotically catenoidal, more precisely it is described as the graph 

2/3 = [a~ 1 a,j + Pj) logr + 0(1) as r — > oo. 

Note that, just as in the minimal surface case, the coefficients of the ends are 
balanced in the sense that they add up to zero. 

It is clear that if two ends are parallel, say flj+i = <ij, we need at least that 
Pj+i—flj > 0, for otherwise the ends of this zero level set would eventually intersect. 
We recall that our further condition on these numbers is that these ends in fact 
diverge at a sufficiently fast rate: 

/3j+i — Pj > 4 maxjal 1 , cr^ 1 } if dj+i = dj ■ (3-4) 

We will explain later the role of this condition. Let us evaluate the error of approx- 
imation S(uq). Using Lemma \2. 2 1 and the fact that w" + f(w) — 0, we find 

S(u ) := A x u + f(u ) = 
-a 2 [\A\ 2 hi + A M hi] w' + 
-a 2 \A\ 2 tw' + 2 a 2 a°- d l h d j h w" + 
a 2 a% (2d t h d j h 1 + d l h 1 d h l ) w" + 
2a 3 (t + h + hx)al j di(h + hx)d j {h + h x )w" + 

a 3 (t + h + h^bld.^ho + w' + a 3 (t + h + h^blw' (3.5) 

where the formula above has been broken into "sizes" , keeping in mind that ho 
is fixed while h x = 0(a). Since we want that uq be as close as possible to be a 
solution of (|f ,ip . then we would like to choose h\ in such a way that the quantity 
(|3.5p be as small as possible. Examining the above expression, it does not look like 
we can do that in absolute terms. However part of the error could be made smaller 
by adjusting hi- Let us consider the "L 2 -projection" onto w'(t) of the error for 
each fixed y, given by 



/oo 
S(u )(y,t) w'(t)dt 
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where for now, and for simplicity we assume the coordinates are defined for all t, 
the difference with the integration is taken in all the actual domain for t produces 
only exponentially small terms in oT x . Then we find 

/OO r-OO 
w' 2 dt + c?di{hQ + h x ) I b\{t + h + h x )w' 2 dt + 
-oo J — oo 

/oo pOO 
(t + h + h)al j w"w'dt + a 3 / (t + h + hi) 3 blw' 2 dt 
-oo J — oo 

(3-6) 

where we have used J_ tw' dt = J_ w"w' dt = to get rid in particular of the 
terms of order a 2 . 

Making all these "projections" equal to zero amounts to a nonlinear differential 
equation for h of the form 

J{h 1 ) = A M h 1 + h 1 \A{y)\ 2 = G {h 1 ) yeM (3.7) 

where Go is easily checked to be a contraction mapping of small constant in hi, in 
the ball radius 0(a) with the C 1 norm defined by the expression in the left hand 
side of inequality (|3.3| . This is where the nondegeneracy assumption on the Jacobi 
operator J enters, since we would like to invert it, in such a way to set up equation 
(|3.7|) as a fixed point problem for a contraction mapping of a ball of the form (|3.3p . 

3.2. Improvement of approximation. The previous considerations are not suf- 
ficient since even after adjusting optimally h, the error in absolute value does not 
necessarily decrease. As we observed, the "large" term in the error, 

-a 2 \A\Hw' + a 2 a° ij d l h Q d 3 h w" 

did not contribute to the projection. In order to eliminate, or reduce the size of 
this remaining part 0(a 2 ) of the error, we improve the approximation through the 
following argument. Let us consider the differential equation 

which has a unique bounded solution with ipo(0) = 0, given explicitly by the formula 
<iPo(t)=w'(t) [ w'(t)- 2 [ sw'(s) 2 ds. 

JO J -oo 

Observe that this function is well defined and it is bounded since J_ sw 1 (s) 2 ds = 
and w'(t) ~ e" <T± l*l as t — > ±oo, with a± > 0. Note also that ipi(t) = \tw'(t) solves 

^) + /'Mt))V>i(t) = «/'(*). 

We consider as a second approximation 

Ux=Uo+(j)\, c/)i(y,t) := a 2 \A(ay)\ 2 ip a (t) ~ a 2 a° j d i h djh (ay)ipi(t) . (3.8) 
Let us observe that 

S(u + <P) - SM + A^ + fM^ + iVo^), TVo(^) = f(u + <t>)-f(u )-f(u )<t>. 
We have that 

dtt&i + f'(uo)<pi = a 2 \A(ay)\ 2 tw' - a 2 a°jdihodjho(ay) w" . 
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Hence we get that the largest remaining term in the error is canceled. Indeed, we 
have 

S( Ul ) = S(u ) - (2a 2 a l3 d l h a d 1 h a w" - a 2 \A(ay)\Hw') + [A x - +JV o (0i). 

Since <p\ has size of order a 2 , a smooth dependence in ay and it is of size 0(r~ 2 e~ ''* ) 
using Lemma [2T2l we readily check that the "error created" 

[A x - Sulfa + MM ■= -a 4 ( l^| 2 #o - a^dihodjho ttf x ) Ah, + R 

satisfies 

\Ro(y,t)\ <Ca 3 (l + r a (y))- 4 e-^. 

Hence we have eliminated the /ii-independent term 0(a 2 ) that did not contribute 
to the projection H(y), and replaced it by one smaller and with faster decay. Let 
us be slightly more explicit for later reference. We have 

S(ui) := A«i + /(«i) = 
-a 2 p| 2 /ii + A M fti] w' + a 2 a°- {d^d^ + di^djho + 8^8^ ) w" 

- a 4 ( \A\ 2 tiP' - a^dihodjho ty{ ) A u h x + 2a 3 {t + h)a\j dihd 3 hw" + Ri (3.9) 
where 

Ri = Ri(y,t, hi(ay), V M hi(ay)) 

with 

|Afii(l/,*,i,j)| + \D 3 Ri(y,t,h3)\ + \Ri(y,t,hj)\ < Ca 3 (l+r a (y))- 4 e~^ 
and the constant C above possibly depends on the number K, of condition (|3 . 3|) . 

The above arguments are in reality the way we will actually solve the problem: 
two separate, but coupled steps are involved: (1) Eliminate the parts of the error 
that do not contribute to the projection n and (2) Adjust h\ so that the projection 
n becomes identically zero. 

3.3. The condition of diverging ends. Let us explain the reason to introduce 
condition (|3 .4(1 in the parameters (3j. To fix ideas, let us assume that we have two 
consecutive planar ends of M, Mj and Mj+i, namely with cij = Oj+i and with 
d = bj+i — bj > 0. Assuming that the normal in Mj points upwards, the coordinate 
t reads approximately as 

t — X3 — oT x bj — h near Mj a , t — a~ 4 bj + \ — 2:3 — h near Mj + i a . 

If we let ho = both on Mj a and Mj + \ a which are separated at distance d/a, then a 
good approximation in the entire region between Mj a and Mj + \ a that matches the 
parts of w(t) coming both from Mj and Mj+i should read near Mj approximately 
as 

w(t) +w(a~ 1 d-t) - 1. 
When computing the error of approximation, we observe that the following addi- 
tional term arises near Mj a : 

E := /( w(t) + w(a~ 1 d — t) — 1 ) - f{w(t)) - f{w(a^d - t)) ~ 

~ [/>(«)) -/'(!)] (wia-'d-t)-!). 
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Now in the computation of the projection of the error this would give rise to 

/oo 
[f'(w(t))-f'(l)](w(a- 1 d~t) - l)w'(t)dt ~c*e-°+i. 
-oo 

where c* ^ is a constant. Thus equation (|3.7p for h\ gets modified with a term 
which even though very tiny, it has no decay as \y\ — > +00 on Mj, unlike the others 
involved in the operator Go in (|3 . T[) . That terms eventually dominates and the 
equation for h\ for very large r would read in Mj as 

A M /i! ~ e-5 ^ 0, 

which is inconsistent with the assumption that h is bounded. Worse yet, its solution 
would be quadratic thus eventually intersecting another end. This nuisance is fixed 
with the introduction of ho satisfying condition (|3.4p . In that case the term E 
created above will now read near Mj a as 

E ~ Ce-°+^ e -(ft+i-ft)iogr„ e - a \t\ = (e-tr-* e -*l*l) 

which is qualitatively of the same type of the other terms involved in the compu- 
tation of the error. 

3.4. The global first approximation. The approximation ui(x) in (|3.2[) will be 
sufficient for our purposes, however it is so far defined only in a region of the type M 
which we have not made precise yet. Since we are assuming that M a is connected, 
the fact that M a is properly embedded implies that R 3 \M Q consists of precisely two 
components S- and S+. Let us use the convention that v points in the direction 
of 5 + . Let us consider the function H defined in R 3 \ M a as 

■w-{-i <«»» 

Then our approximation Ui (x) approaches H(x) at an exponential rate 0(e _<T± ) as 
\t\ increases. The global approximation we will use consists simply of interpolating 
ui with H sufficiently well-inside R 3 \ M a through a cut-off in \t\. In order to 
avoid the problem described in i)3.3l and having the coordinates (y,t) well-defined, 
we consider this cut-off to be supported in a region y-dependent that expands 
logarithmically in r a . Thus we will actually consider a region Afs expanding at 
the ends, thus becoming wider as r a — > 00 than the set Af" previously considered, 
where the coordinates are still well-defined. 

We consider the open set O in M a x M defined as 

O = {(y,t) £M a xR, \t+hx{ay)\ < - + 4 max {a! 1 , ct; 1 } log(l+r Q (y)) =: p a (y) } 

a 

(3.11) 

where 5 is small positive number. We consider the the region AT =: Afs of points x 
of the form 

x = X h (y, t) = y+(t + h (ay) + hx(ay)) v{ay), (y, t) E O, 

namely Ms = Xh{0). The coordinates (y,t) are well-defined in Ms for any suffi- 
ciently small 5: indeed the map Xh is one to one in O thanks to assumption (|3.4j) 
and the fact that hi = 0{a). Moreover, Lemma [2TT1 applies in Ms- 
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Let r](s) be a smooth cut-off function with ri(s) — 1 for s < 1 and = for s > 2. 
and define 

n M — J 7 l(\t + hi(ay)\~ p a (y)~3) if x G Afg , r „ 19 s 

%w - \ o ifx^M 

where p a is defined in (|3.11[) . Then we let our global approximation w(ar) be simply 
defined as 

w := mui + (1 - t? 4 )H (3.13) 
where H is given by (|3. 10[) and u±(x) is just understood to be H(x) outside Afg. 

Since H is an exact solution in R 3 \ Mg, the global error of approximation is 
simply computed as 



5(w) = Aw + /(w) = r) S S(ui) + E (3.14) 

where 

E = 2Vt7 4 Vui + A Vs { Ul - H) + + (1 - %)H) ) - »jj/(ui) . 

The new error terms created are of exponentially small size 0(e~° ) but have in 
addition decay with r a . In fact we have 

\E\<Ce--«r'\ 

Let us observe that \t + hi(ay)\ — \z — ho(ay)\ where z is the normal coordinate 
to M a , hence r]g does not depend on hi, in particular the term Arjg does involves 
second derivatives of h\ on which we have not made assumptions yet. 



4. The proof of Theorem^ 

The proof of Theorem [1] involves various ingredients whose detailed proofs are 
fairly technical. In order to keep the presentation as clear as possible, in this section 
we carry out the proof, skimming it from several (important) steps, which we state 
as lemmas or propositions, with complete proofs postponed for the subsequent 
sections. 

We look for a solution u of the Allen Cahn equation (|1.1| in the form 

u = w + if (4-1) 

where w is the global approximation defined in (|3.13|) and (p is in some suitable 
sense small. Thus we need to solve the following problem 

A(p + f'(v)<p = -S(-w)- N((p) (4.2) 

where 

AT(^) = /( W + ^)-/(w)-/ / (wV 

Next we introduce various norms that we will use to set up a suitable functional 
analytic scheme for solving problem (|4.2|) . For a function g(x) defined in R 3 , 1 < 
p < +oo, /i > 0, and a > we write 

IMLm,* := sup(l + r(ax))> 1 \\g\\ L p {B{x ^ )) , r(x',x 3 ) = \x'\ . 

xgl 3 
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On the other hand, given numbers /i > 0, < a < min{cr+, C-}, p > 3, and 
functions g(y, t) and <f>{y 1 1) defined in M a x R we consider the norms 

\\g\\ p ,^:= sup r a { y ye^(( \f\ p dV a ) . (4.3) 

(y,t)EM a xR \JB((y,t),l) J 

Consistently we set 

llfflloo,/.,* := sup r ct ( 2 /)^e CT|t| \\f\\ L °°(B{{y,t),i)) (4.4) 

(y,t)£M a xl 

and let 

||0||2,p,/i,<r := \\D 2 (j)\\p.jj,.,a + ||-D0||oo,/i,<T + H^HoOj/i.ff ■ (4.5) 

We consider also for a function g(y) defined in M the L p -weighted norm 

WfWnfi ■= [j M \f{y)ni + \yf) p dV(y)y P = \\(l + \yf)f\\ LP{M ) (4.6) 

where p > 1 and /3 > 0. 

We assume in what follows, that for a certain constant /C > and p > 3 we have 
that the parameter function h\(y) satisfies 

IMI* := IMU-(M) + 11(1 +r 2 )I>ML<»(JK) + ll^ 2 ^ill P ,4-f < ACa. (4.7) 

Next we reduce problem (|4.2p to solving one qualitatively similar (equation (|4.20p 
below) for a function <f>(y, t) defined in the whole space M a x K. 

4.1. Step 1: the gluing reduction. We will follow the following procedure. Let 
us consider again ?y(s), a smooth cut-off function with r)(s) = 1 for s < 1 and = 
for s > 2, and define 

Ux) ;= f rtlt + ^oy)!-! +n) if* 6^ _ (48) 



if x g A/i 

We look for a solution of problem (|4. 2(1 of the following form 

<p(x)=&(x)4(y,t)+il>(x) (4.9) 

where <fi is defined in entire M a x R, ^(x) is defined in K 3 and C,2{x)4>(y ,t) is 
understood as zero outside Ns- 

We compute, using that C2 • Ci = Cu 

S(v + <p) = + f'(v)<p + N(tp) + S(v) = 
C 2 [A0 + /'(ui)0 + Ci(/'(«i)+S(t))^ +CiA r W + 0) + ^K)] + 

A^-[(i-Ci)/ , M + Ciff(*)h/' + 

(1 - C 2 )5(w) + (1 - CiM</> + C20) + 2VCiV0 + 0ACi (4.10) 
where H{t) is any smooth, strictly negative function satisfying 

jj (f , _//'(+!) if * > 1, 
^ w ~ 1 /'(-l) if t < -1 • 
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Thus, we will have constructed a solution <p — ^'P + ip to problem (14. 2| if we require 
that the pair (<fi, tji) satisfies the following coupled system 

A0+/ / (ui)0 + Ci(/'(«i)-fr(W+CxW+^) + S(«i)=Ofi)r \t\ < -+3 (4.11) 

a 

AV> + [(i-&)/'(ui) + Off(t)]V' + 

(1 - d)S(w) + (1 - Ci) W + C 2 0) + 2VCiV</» + 0ACi =0 in R 3 . (4.12) 

In order to find a solution to this system we will first extend equation (|4.11[) to 
entire M a x R in the following manner. Let us set 

B(<£) = C 4 [A a - d tt - A a , M J 4> (4.13) 
where A^ is expressed in (y,t) coordinates using expression (|2.14j) and B(0) is 
understood to be zero for \t + hi\ > — + 5. The other terms in equation (|4.11[) are 
simply extended as zero beyond the support of Ci ■ Thus we consider the extension 
of equation (|4.11[) given by 

d a <t> + A„, Mt > +B(0) + f'(w{t))<p = -Sim) 

- - f'W\4> + Cl(/'(«i) - H(t))1> + CiN(i> + </>)} in G M a x R, (4.14) 

where we set, with reference to expression p.9[) , 

§(ui) = -a 2 [L4| 2 /n + A M hi] w' + a 2 a?- {2d l h Q d 3 h 1 + dihidjhi ) to" 

-a 4 ( |A| 2 *V»o- a « dihodjhotyl ) A/n +C 4 [a 3 (t+/i)a^ d i hd j hw"+R 1 (y,t)] (4.15) 
and, we recall 

i?i = Ri(y,t, h\{ay), V M hi(ay)) 

with 

|Ai2i(»,*,*,j)| + P 3 i2i(l/,*,*,j)|+|iii(l/,t,*,j)| <C« 3 (l + r Q (2/))- 4 e- CT l*l. (4.16) 

In summary S'(ui) coincides with S'(iti) if £4 = 1 while outside the support of 
£4, their parts that are not defined for all t are cut-off. 

To solve the resulting system (|4.12p - (|4.14[) , we find first solve equation (|4.12[) in 
i/j for a given <fi a small function in absolute value. Noticing that the potential [ (1 — 
(i)f'(ui) + (iH(t) ] is uniformly negative, so that the linear operator is qualitatively 
like A — 1 and using contraction mapping principle, a solution ip = ^((f>) is found 
according to the following lemma, whose detailed proof we carry out in i j7.1.21 

Lemma 4.1. For all sufficiently small a the following holds. Given <f> with 
||0||2,p,p,<r < 1, there exists a unique solution ip — ^(4>) of problem <\4- such that 

HVlIx := \\D 2 i>\\ p<li< * + \\i/>\\w < Ce-^r. (4.17) 
Besides, VP satisfies the Lipschitz condition 
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- #{<h)\\x < - <hhw ■ (4.18) 

Thus we replace i[> = ^(^O in the first equation (]4. so that by setting 

Ufa) :=B(^) + [/'(«i)-/'(u;)]^ + Ci(/ / («i)-ff(*))*W + Citf(*W)+0), (4.19) 

our problem is reduced to finding a solution <p to the following nonlinear, nonlocal 
problem in M a x M. 

<9 ft + Ay t M a <t> + f{w)4> = -S{u x ) - N(0) in M a x M. (4.20) 

Thus, we concentrate in the remaining of the proof in solving equation (|4.20p . As 
we hinted in ^3.21 we will find a solution of problem (|4.20p by considering two steps: 
(1) "Improving the approximation", roughly solving for <j> that eliminates the part 
of the error that does not contribute to the "projections" j\S(U\) + N(^>)]it/ (t)dt, 
which amounts to a nonlinear problem in (f>, and (2) Adjust hi in such a way that 
the resulting projection is actually zero. Let us set up the scheme for step (1) in a 
precise form. 

4.2. Step 2: Eliminating terms not contributing to projections. Let us 

consider the problem of finding a function (j)(y, t) such that for a certain function 
c(y) defined in M a , we have 

d tt (t> + &y,M a 4> = -Sim) - N(0) + c(y)w'(t) in M a x R, 

(4.21) 

<f)(y,t)w'(t)dt = Q, for all y£M a . 

Solving this problem for <j) amounts to "eliminating the part of the error that does 
not contribute to the projection" in problem (|4.20[) . To justify this phrase let us 
consider the associated linear problem in M a x R 



d tt cp + Ay, Ma <t>+ f'(w(t))cp = g(y,t) + c(y)w'(t), for all {y,t) G M a x M, 

<p(y, t) w'(t) dt = 0, for all y e M a . 

(4.22) 

Assuming that the corresponding operations can be carried out, let us multiply 
the equation by w'(t) and integrate in t for fixed y. We find that 



&v,M a / <f>(y,t)w'dt+ / 4>{y,t)[w"' + f(w)w']dt= / gw' + c(y) 



w 



<y) = - jR T', , (4-23) 



The left hand side of the above identity is zero and then we find that 

J n9(y, t)w'd t 
Jr w ' 2dt 

hence a <j) solving problem (|4.22[) . <f> precisely solves or eliminates the part of g 
which does not contribute to the projections in the equation A</> + f'(w)(j) = g, 
namely the same equation with g replaced by g given by 



g{y,t)=g{y,t)- J % yy ;' w'(t). (4.24) 
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The term c(y) in problem (I4.21| has a similar role, except that we cannot find it 
so explicitly. 

In order to solve problem (|4.2ip we need to devise a theory to solve problem 
(|4.22p where we consider a class of right hand sides g with a qualitative behavior 
similar to that of the error S(u±). As we have seen in (14.151) . typical elements in this 
error are of the type 0((1 + r ct (y))~ M e _CT ' t '), so this is the type of functions g{y, t) 
that we want to consider. This is actually the motivation to introduce the norms 
(|4.3|) . (|4.4[) and (|4.5[) . We will prove that problem (|4.22|) has a unique solution <j> 
which respects the size of g in norm f|4. 3|) up to its second derivatives, namely in 
the norm (|4.5p . The following fact holds. 

Proposition 4.1. Given p > 3 ; (i > and < a < min{er_, cr+}, there exists a 
constant C > such that for all sufficiently small a > the following holds. Given 
f with ||<?||p.^.(T < +oo, then Problem y.g^p with c(y) given by ((713), has 

a unique 

solution (j) with \\(f>\\oo,fi,(T < +oo- This solution satisfies in addition that 

Hh m < C\\g\\ v ^ a . (4.25) 

We will prove this result in §3] . After Proposition 14.11 solving Problem (|4.2ip 
for a small <f> is easy using the small Lipschitz character of the terms involved in the 
operator N(0) in (|4.19p and contraction mapping principle. The error term S(u\) 
satisfies 

\\S{ Ul ) + a 2 Ahiw'\\ pAi(T < Ca 3 . (4.26) 

Using this, and the fact that N(<^) defines a contraction mapping in a ball center 
zero and radius 0(a 3 ) in || ||2,p,4,o-, we conclude the existence of a unique small 
solution <j) to problem (|4.2ip whose size is 0(a 3 ) for this norm. This solution <fi 
turns out to define an operator in hi <ft = ^(/ii) which is Lipschitz in the norms || ||« 
appearing in condition (|4.7p . In precise terms, we have the validity of the following 
result, whose detailed proof we postpone for §7.21 

Proposition 4.2. Assume p>3,0</i<3, < a < min{a+,cr_}. There exists 
a K > such that problem |7.&P has a unique solution <f> — ^(/ii) such that 

Besides, $ has a Lipschitz dependence on hi satisfying l\4- 7| in the sense that 

||*(/ii) - 3>(h 2 )h, P ,w < Ca 2 \\hi - h 2 \\*. (4.27) 

4.3. Step 3: Adjusting hi to make the projection zero. In order to conclude 
the proof of the theorem, we have to carry out the second step, namely adjusting hi, 
within a region of the form (|4.7p for suitable JC in such a way that the "projections" 
are identically zero, namely making zero the function c(y) found for the solution 
<j) = Q(hi) of problem (I4.2ip . Using expression (|4.23p for c(y) we find that 

c( y ) I w ' 2 = [ S(ui)w'dt+ [ N($(/ii))u/di. (4.28) 
Jr Jr Jr 

Now, setting c* := J R w' 2 dt and using same computation employed to derive 
formula (|3.6p . we find from expression (|4.15[) that 



/ S(u 1 )(y,t)w'(t)dt = -c*a 2 {& M hi + hi\A\ 2 ) + c*a 2 G*i(/ii) 
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where 

c*Gi(fti) = -a 2 A/ii(|A| 2 / tif)' w' dt - a^dihodjho / ttp^w' dt) + 

Jr Jwl 

adi(ho + hi)dj(ho + hi) / ^4,(t + h)a i jw"w' dt +a~ 2 / (4 Ri(y, t, h\, V Aih\ ) w' dt 
Js. Jm 

(4.29) 

and we recall that R\ is of size 0(a 3 ) in the sense (|4.16p . Thus, setting 

c,G 2 (/ii) := a' 2 [ N($(hi))w'dt, G{h x ) := Gi{hx) + G 2 (h 1 ), (4.30) 



we find that the equation c(y) = is equivalent to the problem 

J(hx) = A M hi + \A\ 2 h x = G(/n) in M. (4.31) 

Therefore, we will have proven Theorem [1] if we find a function h\ defined on 
M satisfying constraint (|4.7[) for a suitable JC that solves equation (|4.31|) . Again, 
this is not so direct since the operator J has a nontrivial bounded kernel. Rather 
than solving directly (14.3 1 [) . we consider first a projected version of this problem, 
namely that of finding h\ such that for certain scalars C\,...,cj we have 

,/ 

J(h 1 )=G(h 1 )+Y'-^Zi inM, 



M 1 



^—^dV = 0, i = l,...J. (4.32) 



Here Zi,...,zj is a basis of the vector space of bounded Jacobi fields. 

In order to solve problem (|4. 32|) we need a corresponding linear invertibility 
theory. This leads us to consider the linear problem 

J 



J(h) = f + J2rf^ 2 * inA/ ' 
i=i 

dV = 0, i = l,...J. (4.33) 



i=l 

Zih 



M 1 



Here z\, ...,zj are bounded, linearly independent Jacobi fields, and J is the dimen- 
sion of the vector space of bounded Jacobi fields. 



We will prove in §6. II the following result. 

Proposition 4.3. Given p > 2 and f with H/IL44 < +00, there exists a unique 
bounded solution h of problem |^T55| . Moreover, there exists a positive number 
C = C{p, M) such that 

\\h\U := ||ft,||oo + II (1 + \y?) DhU + || D 2 h \\ pA _ ± < C||/|| p!4 _4 . (4.34) 

Using the fact that G is a small operator of size 0(a) uniformly on functions 
hi satisfying (|4.7p , Proposition 14.31 and contraction mapping principle yield the 
following result, whose detailed proof we carry out in ^JHJ 

Proposition 4.4. Given p > 3, there exists a number K, > such that for all 
sufficiently small a > there is a unique solution h\ of problem 32\ that satisfies 
constraint 
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4.4. Step 3: Conclusion. At the last step we prove that the constants Ci found 
in equation (14.32p are in reality all zero, without the need of adjusting any further 
parameters but rather as a consequence of the natural invariances of the of the full 
equation. The key point is to realize what equation has been solved so far. 

First we observe the following. For each h\ satysfying (|4.7p . the pair (0, tjj) with 
<j> = Q(hi), ijj = ^f(4>), solves the system 

A<j>+f'(ui)<f> + Ci(/'(ui)-#W)V +ClN(ip + <t>) + S( Ul ) = c(y)w'(t) for \t\ < - + 3 

a 

A^ + [(i-Ci)/'K) + Ciff(*)]V> + 

(1 - Cj)5(w) + (1 - (i)N(ip + C, 2 <t>) + 2VC1V0 + 0ACi - in R 3 . 
Thus setting 

<f(x) = (2(x)(j>(y,t) + ip(x), u = w + if, 
we find from formula (|4.10j) that 

Au + f(u) = S(w + <p) = ( 2 c(y) w'(t) . 
On the other hand choosing hi as that given in Proposition ^. 4l which solves problem 
(|4.32p . amounts precisely to making 

for certain scalars Cj. In summary, we have found h\ satisfying constraint (|4.7|) such 
that 

u = w + 1 (a;)$(ft 1 ) + *($(/n)) (4.35) 

solves the equation 

J 

Au + f(u)=Y,y^-rUayW(t) (4.36) 

where Sj = c*a 2 Ci. Testing equation (|4.36|) against the generators of the rigid 
motions diU i = 1,2,3, —X2d\U + Xi&2U, and using the balancing formula for the 
minimal surface and the zero average of the numbers (3j in the definition of ho, we 
find a system of equations that leads us to a = for all i, thus conclude the proof. 
We will carry out the details in fj9l 

In sections ^SJHl we will complete the proofs of the intermediate steps of the 
program designed in this section. 

5. The linearized operator 

In this section we will prove Proposition 14.11 At the core of the proof of the 
stated a priori estimates is the fact that the one-variable solution w of is 
nondegenerate in L°°(1R 3 ) in the sense that the linearized operator 

L{(j>) = A y( j) + du(t) + f'(w(t))cf>, {y, ()eI 3 =R 2 x M, 

is such that the following property holds. 
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Lemma 5.1. Let cf> be a bounded, smooth solution of the problem 

L{cj)) = in R 2 x K. (5.1) 
Then <j>(y,t) = Cw'{t) for some Cel. 

Proof. We begin by reviewing some known facts about the one-dimensional operator 
Lo{4>) = ip" + f'(w)ip. Assuming that tp(t) and its derivative decay sufficiently fast 
as \t\ — > +oo and defining ip{t) = w'(t)p(t), we get that 

-f'(w)^ 2 ]dt = { L (i>)i>dt = fw' 2 \p'\ 2 dt, 



therefore this quadratic form is positive unless ip is a constant multiple of w' . Using 
this and a standard compactness argument we get that there is a constant 7 > 
such that whenever J R tpw' — with ijj S we have that 

( |V'] 2 - /'M^ 2 )dt> 1 j{ |Vf + H 2 ) dt. (5.2) 



Now, let 4> be a bounded solution of equation (15. ip . We claim that tfi has exponential 
decay in i, uniform in y. Let us consider a small number a > so that for a certain 
to > and all |t| > i we have that 

f'(w) < -2a 2 . 

Let us consider for e > the function 

2 

g e (t, y) = e — (l*l-*o) + £ £ cosh^) 
Then for |t| > to we get that 

L(^) < if |i| > t . 
As a conclusion, using maximum principle, we get 

|0| < U\\oo9e if|t|>t , 

and letting e — ► we then get 

Wtf,t)| < Cll^lUe^l*! if |t| >t . 
Let us observe the following fact: the function 

">'(*) 



<j>(y,t) = <j>(y,t) - U u/(C)#y,CK 
also satisfies L(0) = and, in addition 



w'(t)4>(y,t)dt = for all y e R 2 . (5.3) 
In view of the above discussion, it turns out that the function 

<p{y) ■= [ 4> 2 iv,t)dt 



is well defined. In fact so are its first and second derivatives by elliptic regularity of 
</>, and differentiation under the integral sign is thus justified. Now, let us observe 
that 



A y <p(y) = 2 A y 4> -4>dt + 2 / \V y t 

JR JR 
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and hence 

o = / • 4>) 



(5.4) 

= -A yip - / \v y 4>\ 2 dz- / (\4> t \ 2 -f'H4> 2 )dt. 

1 Jr Jr 

Let us observe that because of relations (15.31) and (15.21) , we have that 



(i<H 2 -/V)<ncft >if. 



It follows then that 

1 

2" 

Since tp is bounded, from maximum principle we find that ip must be identically 
equal to zero. But this means 



Aylfi — Jlfi > 0. 



4>{y, t) = (f w'(C) 4>{y, C) dc] (5.5) 



Then the bounded function 

g(y)= I ™c(0Hy,t)d( 

Jr 

satisfies the equation 

A y g = 0, in K 2 . (5.6) 

Liouville's theorem implies that g = constant and relation (15. 5p yields <fi(y, t) = 
Cw'(t) for some C. This concludes the proof. □ 

5.1. A priori estimates. We shall consider problem (|4.22|) in a slightly more 
general form, also in a domain finite in y-direction. For a large number R > let 
us set 

M* := {y G M a / r(ay) < R} 
and consider the variation of Problem ()4.22j) given by 



(5.7) 



d tt <j> + Ayj\j a tfi + f'{w{t))dp = g{y, t) + c(y)w'(t) in x R, 

0=0, on dM* x E, 

<f>(y, t) w'(t) dt = for all y G M*, 

where we allow R = +oo and 

c(y) / w' 2 dt = - / g(y,t)w'dt. 
Jr Jr 

We begin by proving a priori estimates. 

Lemma 5.2. Let us assume that < a < min{ <r_, a+} and \i > 0. Then there 
exists a constant C > such that for all small a and all large R, and every solution 
<j) to Problem (|5. J3[) with \\4>\\oo.^,<y < +oo and right hand side g satisfying \\g\\ p ,fj,,a < 
+oo we have 

ll-D 2 0IW + \\D<t>\\ < C\\g\\ Pill „„. (5-8) 
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Proof. For the purpose of the a priori estimate, it clearly suffices to consider the 
case c(y) = 0. By local elliptic estimates, it is enough to show that 

HWoo^a < C\\g\\ p .^. (5.9) 

Let us assume by contradiction that (|5.9[) does not hold. Then we have sequences 
a = a n —> 0, R = R n —> oo, g n with \\g n \\p^,cr —> 0, <f> n with \\<p n \\ = 1 such 
that 



d tt (pn + &y,M a <t>n + f (w{t))<j> n = g n in X R, 

d> n = on dM* x R, 



(5.10) 



4> n (y, t) w' (t) dt = Q for all y £ M* . 
Then we can find points (y n ,t n ) S x R such that 

e-°^\l+r{a n y n )y\<l> n {y n ,t n )\ > i. 

We will consider different possibilities. We may assume that either r a {y n ) — 0{\) 
or r a (y n ) -> +oo. 

5.1.1. Case r(a n y n ) bounded. We have a n y n lies within a bounded subregion of M, 
so we may assume that 

a n y n -> yo e M. 

Assume that j/o € Yk(Uk) for one of the local parametrization of M. We consider 
y„, y € % with Y k (y n ) = a n y n , Y k (y ) = y . 

On a~ x Y k {U k ), M a is parameterized by Yk, a „(y) = a^ 1 y fe (a„y), y G a^Uk- 
Let us consider the local change of variable, 

y = o^yn + y- 

5.1.2. Subcase t n bounded. Let us assume first that \t n \ < C. Then, setting 

4>n(y, t) := 4> n (a~ 1 y n + y, t), 

the local equation becomes 

a%{y n + a n y)9ij4> n + antfifn + a n y)dj4> n + d tt 4> n + f'(w(t))4> n = g n {y,t) 

where g n (y, t) := g n (yn + cty, t). We observe that this expression is valid for y well- 
inside the domain a~ l Uk which is expanding to entire R 2 . Since (f) n is bounded, and 
g n — ► in Lf oc (M. 2 ), we obtain local uniform iy 2 ' p -bound. Hence we may assume, 
passing to a subsequence, that <j> n converges uniformly in compact subsets of M 3 to 
a function <p(y, t) that satisfies 

a?j(y)M + ^ + /'K*))^ = 0. 

Thus 4> is non-zero and bounded. After a rotation and stretching of coordinates, 
the constant coefficient operator a^(y)dij becomes A y . Hence Lemma \b . 1 1 implies 

that, necessarily, <f>(y,t) — Cw'(t). On the other hand, we have 

4> n (y, t) w' (t) dt — > / 4>(y, t) w' (t) dt asn— >oo. 
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Hence, necessarily <j> = 0. But we have (1 + r(a„y„)) M \4> n (0,t n )\ > |, and since t n 
and r(a n y n ) were bounded, the local uniform convergence implies 0^0. We have 
reached a contradiction. 

5.1.3. Subcase t n unbounded. If y n is in the same range as above, but, say, t n — > 
+oo, the situation is similar. The variation is that we define now 

K (y, *) = e CT(t " +t) <£„ (a" 1 y„ + y , t n + t) , g n (y, t) = e a (t " +t) .g„ (a^Vn + y, + 1) . 
Then n is uniformly bounded, and g n — > in L]' oc (]R 3 ). Now </> n satisfies 
0°j(yn + "«y) %0n + 9 t t4>n + ctnbjijn + a n y) dj(j) n 

-2adt4> n + (f'(w(t + t n ) + cr 2 )4> n = g n . 
We fall into the limiting situation 

a*j dij4> + du4> ~ 2cr d t 4> - {a\ - a 2 ) = in M 3 (5.11) 

where a*j is a positive definite, constant matrix and 4> ^ 0. But since, by hypothesis 
crL — a 2 > 0, maximum principle implies that cj) = 0. We obtain a contradiction. 

5.1.4. Case r(a n y n ) — > +oo. In this case we may assume that the sequence a n y n 
diverges along one of the ends, say Mk- Considering now the parametrization 
associated to the end, y — ipk(y), given by (|2.ip . which inherits that for M ani fc, 
y = a~ 1, )/'fc(a„y). Thus in this case (y„ + a n y) — > <5y, uniformly in compact 
subsets of M 2 . 

5.1.5. Subcase t n bounded. Let us assume first that the sequence t n is bounded and 
set 

4>n(y,t) = (1 + r(y n + a n y)Y (j) n (a~ l y n + y,t n + t). 

Then 

dj(r~^4> n ) = -nar~ tJ, ~ 1 djr<t> + r^^djtp 
dij{r-y n ) = fi(ji + l)a 2 r-^ 2 d l rd J r^ - /ioV^fiyr^ - m^^djrd^ 

Now diT = 0(1), dijT = 0(r _1 ), hence we have 

dj(r-y n )=r-» [d^ + 0(ar~ 1 ), 

9ij(r-y n ) = r"" + 0(ar~ 1 )d^ + 0(a 2 r~ 2 ) 

and the equation satisfied by <f> n has therefore the form 

A y </>„ + d t t4>n + o(l)dij4> n + o(l)dj(f> n + o(l) 4> n + f(w(t))4> 

where </>„ is bounded, g n — ► in Lf oc (R 3 ). From elliptic estimates, we also get 
uniform bounds for ||<9j</>„||oo and ||9ij0n||p,o,o- In the limit we obtain a 4> ^ 
bounded, solution of 

A y 4> + d tt + f'(w(t))4> = O, f $(y,t)w'(t)dt = 0, (5.12) 

a situation which is discarded in the same way as before if <f> is defined in R 3 . There 
is however, one more possibility which is that r(a n y n ) — R n = 0(1). In such a case 
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we would see in the limit equation (|5.12j) satisfied in a half-space, which after a 
rotation in the y-plane can be assumed to be 

H = {(y,t) e R 2 x R / y 2 < 0}, with 0(y x , 0, t) = for all (yi,t)el 2 . 

By Schwarz's reflection, the odd extension of <f>, which achieves for 7/2 > 0, 
<fi(jji,y2,t) = —(f>(jji 1 —y2,t) 1 satisfies the same equation, and thus we fall into one 
of the previous cases, again finding a contradiction. 

5.1.6. Subcase t n unbounded. Let us assume now \t n \ — > +00. If t n — > +00 we 
define 

4>„(y,t) = (1 + r(y„ + a B y))" e*" +t M^fn + y, t„ + t). 
In this case we end up in the limit with a </> 7^ bounded and satisfying the equation 
A y 4> + d tt 4> - 2ad t 4> - (<4-ct 2 )0 = 

cither in entire space or in a Half-space under zero boundary condition. This implies 
again cf> — 0, and a contradiction has been reached that finishes the proof of the a 
priori estimates. □ 

5.2. Existence: conclusion of proof of Proposition 14.11 Let us prove now 
existence. We assume first that g has compact support in M a x R. 

d tt + \,M«<t> + f\w{t))(t) = g(y, t) + c(y)w'(t) in x R, 

0=0, on dM* x R, 

(f)(y, t) w'(t) dt = for all y € M^, 



(5.13) 



-00 

where we allow R = +00 and 

c(y) / w' 2 dt = — / g(y,t)w'dt. 
Jr Jm 

Problem (|5.13[) has a weak formulation which is the following. Let 

H = {(j> e Hl{M* X M) / [ 4>(y,t)w'(t)dt = for all y e M*} . 



H is a closed subspace of Hq(M^ x R), hence a Hilbert space when endowed with 
its natural norm, 



2 

\H 



( |5 t 0| 2 + |V M >| 2 - /'M*) 2 ) rfF a * 
is then a weak solution of Problem (|5.13|) if (/) G H and satisfies 

o(0,V) := / (V M «f Vm^ - /'(«?(*)) 0V) dKxdt 

JM»xR 



g V rfV'a d< for all ?/> e H. 

It is standard to check that a weak solution of problem (|5.13[) is also classical 
provided that g is regular enough. Let us observe that because of the orthogonality 
condition defining H we have that 

7/ %l?dV a dt < a(V>,V0 for all if; € H. 
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Hence the bilinear form a is coercive in H, and existence of a unique weak solution 
follows from Riesz's theorem. If g is regular and compactly supported, tp is also 
regular. Local elliptic regularity implies in particular that <p is bounded. Since for 
some to > 0, the equation satisfied by <fi is 

+ f'{w(i)) <t> = c(y)w'(t), \t\>t , y€M*, (5.14) 

and c(y) is bounded, then enlarging to if necessary, we see that for a < min{(7_|_, er-}, 
the function v(y,t) := Ce^ a ^ + ee cr ' t is a positive supersolution of equation (|5. 14[) . 
for a large enough choice of C and arbitrary e > 0. Hence |0| < Ce~ <T ' t ', from 
maximum principle. Since is bounded, we conclude that H^Hp^.o- < +oo. 
From Lemma [5721 we obtain that if R is large enough then 

00,/J-,CT H~ || ^|| 00,/i.,£T 

< C\\g\\ p ^ a (5.15) 

Now let us consider Problem (|5.13|) for R = +oo, allowed above, and for ||<7|| P) ^ )CT < 
+oo. Then solving the equation for finite R and suitable compactly supported 
gn, we generate a sequence of approximations 4>R which is uniformly controlled in 
R by the above estimate. If g^ is chosen so that gn — > g in L^ oc (M a x R) and 
||<?i?||p,/i,<T < C||<?||p,/i,<T, We obtain that 4>r is locally uniformly bounded, and by 
extracting a subsequence, it converges uniformly locally over compacts to a solution 
<j) to the full problem which respects the estimate (|4.25j) . This concludes the proof 
of existence, and hence that of the proposition. □ 

6. The Jacobi operator 

We consider this section the problem of finding a function h such that for certain 
constants c\, . . . ,cj, 

J 

J(h) = A M h+\A\ 2 h = f + Y^-^z t inM, (6.1) 

and prove the result of Proposition [43] We will also deduce the existence of Jacobi 
fields of logarithmic growth as in Lemma I3TT1 We recall the definition of the norms 

II \\ P ,p m flU). 

Outside of a ball of sufficiently large radius Rq, it is natural to parameterize each 
end of M, 7/3 = Fk(l/l> 2/2) using the Euclidean coordinates y = (y%, 1/2) € R 2 . The 
requirement in / on each end amounts to / € L p (B(0, 1/Rq)) where 

/(y):=|yr 4 /(|yr 2 y). (6.3) 

Indeed, observe that 

ll/HWoi/*o))= / |yr 4p l/(|yr 2 y)l^y = / |y| 4(p - 1) l/(y)l p dy • 

In order to prove the proposition we need some a priori estimates. 

Lemma 6.1. Let p > 2. For each Rq > sufficiently large there exists a constant 
C > such that if 

Wf\\pA-§ + < +00 
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and h solves 

A M h+\A\ 2 h = f, yeM, \y\>R , 

then 

\\h\\L°°(\ y \>2R ) + II \y\ 2 Dh\\ L ^ (lyl>2Ro) + || \y\ 4 ~pD 2 h\\ LPilyl>2Ro) < 

C[\\f\\p,4,-± + IMU°°(Bo<|2/|<3B ) ]• 

Proof. Along each end Mk of M, Am can be expanded in the coordinate y as 

A M = A + 0(|y|- 2 )£> 2 + 0(|y|- 3 ) J D. 
A solution of h of equation (|6.ip satisfies 

A M h+\A\ 2 h = f, \y\>R 
for a sufficiently large Rq. Let us consider a Kelvin's transform 

My) = %/|y| 2 )- 

Then we get 

A%)^|yr 4 (A/ l )(y/|y| 2 ). 

Besides 

0(\y\- 2 )D 2 h(y) + 0(\y\- 3 )Dh(y) = 0(\y\- 6 )D 2 h(y/\y\ 2 ) + 0(|y|- 5 )£>%/|y| 2 ) . 
Hence 

(A M ^)(y/|y| 2 ) = |y| 4 [ Ah(y) + 0(\y\ 2 )D 2 h(y) + 0(\y\)D~h(y) . 

Then h satisfies the equation 

Ah + 0(\y\ 2 )D 2 h + 0(\y\)Dh + 0(l)h = /(y), < |y| < i- 

Kq 

where / is given by (|6.3[) . The operator above satisfies maximum principle in 
B(0, -g-) if i? is fixed large enough. This, the fact that h is bounded, and L p - 
elliptic regularity for p > 2 in two dimensional space imply that 

\\Hl^(B(0,1/2R )) + H-E > ' l IU° (.B(0,l/2ilo)) + II L»(B(0,l/2B )) - 

C[II/IUp((B(0,1/-Ro)) + IWU°°(l/3Bo<|y|<l/B )] ^ 

c[||/iu_| + \\h\\ L~(B(fl <|y|<3B )) J ■ 

Let us observe that 

l|/l|U°°(B(0,l/2fl )) = INU<*>(|y|>2fl ): 
||-Dh|U»(B(Q,l/2flo)) = II |y| 2 -D^I|L~(|y|>2B )- 

Since 

l^ 2 %)l < C{ |y|- 4 \D 2 h( |y|- 2 y)| + |y|- 3 |Dji( |y| ) 

then 

|y| 4 -^ 2 %)l < C-dyl-^l^dyl-^l + lyl-F^I^MIyl-V)!). 

Hence 



/ |y| 4 f- 4 |Z? 2 ^rdy < 

■%|>2Bo 
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C( I \D 2 h(yWdy + \\Dh\\l^ B{01/2Ro)) [ M^dy). 

Jb(Q,1/2R ) ' J\j\>2R 

It follows that 

l|/«IU-(|y|>2fi ) + II |y| 2 -D^I|L~(|y|>2flo) + II \j f~ ? D 2 h\\ Lv ( | y | >2R()) < 
C[ll/llp,4-| + ll /l llL~(S(fl <|y|<3fl.o)) ] • 

Since this estimate holds at each end, the result of the lemma follows, after possibly 
changing slightly the value Rq. □ 



Lemma 6.2. Under the conditions of Lemma \ 6.1l assume that h is a bounded 
solution of Problem (I6'.il) - (|6'.^l) . Then the a priori estimate 34\ ) holds. 



Proof. Let us observe that this a priori estimate in Lemma 16.11 implies in particular 
that the Jacobi fields z% satisfy 

Vli(y) = 0{\y\~ 2 ) as \y\ -> +co. 
Using Zi as a test function in a ball B(0, p) in M we obtain 

(hd v Zi - Zid v Zi) + / (A M z, + \A\ 2 Zi) h = 

.'/»■((),„> J\v\<p 

J 



f A • » / 

J\v\<o J A 



\y\<p ~ JM^ + r 4 

Since the boundary integral in the above identity is of size 0(p^ 1 ) we get 

/* + £>/ (6.4) 
so that in particular 

\c 3 \<C\\f\\ pA _ ± for all j = l,...,J. (6.5) 

In order to prove the desired estimate, we assume by contradiction that there 
are sequences h n , f n with ||/i n ||oo = 1 and ||/n|L 4_i — > 0, such that 

j n , 

A M h n + \A\ 2 h n = f n + y i T i r 1 J 
3=1 

h n Zi 



= for all i = 1 , . . . , J. 

M 



1 + r 4 



Thus according estimate (16.51) . we have that c" — > 0. From Lemma ROl we find 

ll^n|U~(|a|>2fl ) ^ + \\hn\\L°°(.B(0,3Ro))] ■ 

The latter inequality implies that 

HMU°°(B(0,3ilo)) > 7 > 0- 

Local elliptic estimates imply a C 1 bound for on bounded sets. This implies 
the presence of a subsequence h n which we denote the same way such that h n —* h 
uniformly on compact subsets of M, where h satisfies 

A M h+\A\ 2 h = 0. 
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h is bounded hence, by the nondegeneracy assumption, it is a linear combination 
of the functions Besides h ^ and satisfies 

- ^ , — for all i = 1, . . . , J . 

l + r 4 

The latter relations imply h = 0, hence a contradiction that proves the validity of 
the a priori estimate. □ 



/ 

Jm 



6.1. Proof of Proposition 14731 Thanks to Lemma [6.21 it only remains to prove 
existence of a bounded solution to problem (|6.ip - (|6.2|) . Let / be as in the statement 
of the proposition. Let us consider the Hilbert space H of functions h 6 Hl oc {M) 
with 

\\h\\%:= j \Vhf + -^\h\ 2 <+^, 
f 1 

/ jhh = for all i = 1, .... J. 

Jm 1 + 7,4 

Problem (|6.1|) - (16.2p can be formulated in weak form as that of finding h € H with 

/ VWV> - |v4| 2 /i V = - / fi> for all ^eff. 
^a/ Jm 

In fact, a weak solution h € i? of this problem must be bounded thanks to elliptic 
regularity, with the use of Kelvin's transform in each end for the control at infinity. 
Using that \A\ 2 < Cr~ 4 , Riesz representation theorem and the fact that H is 
compactly embedded in L 2 ((l + r i )~ 1 dV) (which follows for instance by inversion 
at each end), we see that this weak problem can be written as an equation of the 
form 

h - T(h) = f 

where T is a compact operator in H and / € H depends linearly on /. When / = 0, 
the a priori estimates found yield that necessarily h = 0. Existence of a solution 
then follows from Fredholm's alternative. The proof is complete. □ 



6.2. Jacobi fields of logarithmic growth. The proof of Lemma l3.1l We will 
use the theory developed above to construct Jacobi fields with logarithmic growth 
as r — * +oo, whose existence we stated and use to set up the initial approximation 
in Lemma l3Tl One of these Jacobi fields is the generator of dilations of the surface, 
z o(y) = y' v (y)- We will prove next that there are another m— 2 linearly independent 
logarithmically growing Jacobi fields. 

Let us consider an m-tuple of numbers 0i, ■ . ■ ,/3 m with J2j 0j = 0, and any 
smooth function p(y) in M such that on each end Mj we have that for sufficiently 
large r = r(y), 

p(y) = (-iy/3 j logr(y), y £ M 3 

for certain numbers , (3 m that we will choose later. To prove the result of 

Lemma 13.11 we need to find a solution ho of the equation ^J(ho) = of the form 
ho = p + h where h is bounded. This amounts to solving 



J{h) = -J( P ) 



(6.6) 
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Let us consider the cylinder Cr = {x el 3 / r(x) < R} for a large R. Then 
/ J( P ) z 3 dV = [ J{z 3 )z 3 dV + f (z 3 d nP - pd n z 3 ) da(y) . 

JMHCr JMnCn JdCnHM 



Hence 



/ J{p) z 3 dV = V / {z 3 d n p - pd n z 3 ) da(y) . 

JMnC R - =1 JdCnHMj 

Thus using the graph coordinates on each end, we find 

/ J(p)z 3 dV = 

JMnCn 



£(-!) j 

J'=l 



-£ [ v 3 da(y) - Pj log R [ d r v 3 da(y) 

R J\y\=R J\y\=R 



+ OiR- 1 ). 



We have that, on each end M, 



Hence we get 



(-l) j 



Vl+|VF fc (y)|2 



= (-iy +0(r~ 2 ), d r is 3 (y) = 0(r- 3 ). 



/ J( P )z 3 dV = 2TrJ2f3 3 +0(R- 1 ) 

JMnCn =1 



It is easy to see, using the graph coordinates that J{p) = 0(r 4 ) and it is hence 
integrable. We pass to the limit R — > +00 and get 



« m 

/ J( P )z 3 dV = 2ttV/?, = . 



(6.7) 



We make a similar integration for the remaining bounded Jacobi fields. For 
2» = Vi{y) i = 1, 2 we find 

/ J(p)^2rfV = V(-f) J v 2 da( V ) - fologR [ d r v 2 da(y) +0(R- 1 ). 

JMnCn „•_-, J|v|=fl J|vl=« 



Now, on Mj , 



Ml) 



Hence 



/ J(P). 

J M 



idV = i = 1,2. 



Finally, for z 4 (y) = (-t/2, yi, 0) • we find on Mj, 
(-l) j z 4 (y) = -y2d 2 F j +y 1 d 1 F j = b J-^ - b 2 J± + 0( r - 2 ), d r z 4 = 0( r - 2 ) 



and hence again 



/ J(p)z 4 dV = 
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From the solvability theory developed, we can then find a bounded solution to the 
problem 

./ 

3=1 

Since J M J(p)zidV = and hence J M J(p)zidV — 0, relations (|6.4p imply that 
Cj = for all i. 

We have thus found a bounded solution to equation (|6.6p and the proof is con- 
cluded. □ 

Remark 6.1. Observe that, in particular, the explicit Jacobi field zo(y) = y ■ v(y) 
satisfies that 

z(y) = (-l) j aj logr + 0(1) for all y € Mj 

and we have indeed ^ • ctj — 0. Besides this one, we thus have the presence of 
another m — 2 linearly independent Jacobi fields with \z(y)\ ~ logr as r — > +oo, 
where m is the number of ends. 

These are in reality all Jacobi fields with exact logarithmic growth. In fact if 
J(z) = and 

|*(y)|<Clogr, (6.8) 

then the argument in the proof of Lemma |6 . 1 1 shows that the Kelvin's inversion z(y) 
as in the proof of Lemma 16.21 satisfies near the origin Az — f where / belongs to 
any L p near the origin, so it must equal a multiple of log |y| plus a regular function. 
It follows that on My there is a number j3j with 

2(y) = (-imiog|y| + / l 

where h is smooth and bounded. The computations above force J^j Pj = 0- ^ 
follows from Lemma 13.11 that then z must be equal to one of the elements there 
predicted plus a bounded Jacobi field. We conclude in particular that the dimension 
of the space of Jacobi fields satisfying (|6.8[) must be at most m— 1+ J, thus recovering 
a fact stated in Lemma 5.2 of [33]. 



7. Reducing the gluing system and solving the projected problem 

In this section we prove Lemma 14. 1[ which reduces the gluing system (|4.12[1 - 
(|4.14p to solving the nonlocal equation (|4.20|) and prove Proposition ^. 21 on solving 
the nonlinear projected problem (|4.2ip . in which the basic element is linear the- 
ory stated in Proposition 14.11 In what follows we refer to notation and objects 
introduced in fl4~2l 

7.1. Reducing the gluing system. Let us consider equation (|4.12p in the gluing 
system ([d2]) - ([ST3]) . 

AV-W^ Q (a;)^ + (l-C2)^(w) + (l-Ci)A r (^ + C20) + 2VCiV(/) + (/)ACi -0 in M 3 

(7.1) 

where 

W a (x) := [(l-Ci)/'(«i) + &#(*)] ■ 
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7.1.1. Solving the linear outer problem. We consider first the linear problem 

A^> - W a (x)ij) + g(x) = in M 3 (7.2) 

We observe that globally we have < a < W a (x) < b for certain constants a and 
b. In fact we can take a = min{cr^, o+} — r for arbitrarily small r > 0. 

We consider for the purpose the norms for 1 < p < +oo, 

\\g\\p,n ■= sup (1 + KaaO^IMI^BOM)), r ( x '> x a) = \ x '\ ■ 

Lemma 7.1. Given p > 3, /i > 0, there is a C > such that for all sufficiently 
small a and any g with \\g\\ p ^ < +oo there exists a unique ip solution to Problem 
Q7.if with HV'lloo./j < +oo. This solution satisfies in addition, 

P 2 VIU + IMIoc, m < C\\g\\ p ^. (7.3) 

Proof. We claim that the a priori estimate 

||VIU„ < C\\g\\M (7.4) 

holds for solutions if) with ||^||oo,/i < +°° to problem (|7.2p with HgUp,^ < +oo 
provided that a is small enough. This and local elliptic estimates in turn implies 
the validity of (|7.3p . To see this, let us assume the opposite, namely the existence 
a n —>■ 0, and solutions ip n to equation (|7.2|) with HV'ralloo,^ = 1, HflVillp,^ —> 0. Let 
us consider a point x n with 

(1 + r{a n x n )Ytp n {x n ) > ^ 

and define 

tpn(x) = (1 + r(a n (x n + x)) tJ "ip n (x n + x), g n (x) = (1 + r(a n (x n + x)Yg n {x n + x), 

W n (x) = W an {x n +x). 
Then, similarly to what was done in the previous section, we check that the equation 
satisfied by ip n has the form 

Aljj n - W n {x)i> n + o(l)V^ n + o(l)V>n = 9n- 

ip n is uniformly bounded. Then elliptic estimates imply i°°-bounds for the gradient 
and the existence of a subsequence uniformly convergent over compact subsets of 
M. 3 to a bounded solution ip ^ to an equation of the form 

Aijj - W, (x)tp = in W 3 

where < a < W* (x) < b. But maximum principle makes this situation impossible, 
hence estimate (|7.4I) holds. 

Now, for existence, let us consider g with Hgllp,^ < +oo and a collection of 

approximations g n to g with ||g n ||oo, AI < +oo, g n — > g in Lf oc (R 3 ) and ||flri||p,M < 
C||<7|| P)(tl . The problem 

Aifj n - W n (x)?p n = g n inM 3 

can be solved since this equation has a positive supersolution of the form 
C n (l + r(ax) ) _AI , provided that a is sufficiently small, but independently of n. Let 
us call ip n the solution thus found, which satisfies || V , n||oo,^* < +oo. The a priori 
estimate shows that 

PVnIU + HiMU,, < C\\g\\ Ptfi . 
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and passing to the local uniform limit up to a subsequence, we get a solution ip to 
problem (17. 2|) . with ||^||oo,/j < +00. The proof is complete. □ 



7.1.2. The proof of Lemma \4--1\ Let us call ip :— T(a) the solution of Problem (|7.2|) 
predicted by Lemma 17. II Let us write Problem (|7.ip as fixed point problem in the 
space X of W^-functions -0 with \\ip\\x < +00, 

i ! = T(g 1 +K^)) (7.5) 

where 

9i = (1 - C 2 )5(w) + 2VCiV0 + 0ACi , Kty) = (1 - Ci)JV(V> + C20) ■ 
Let us consider a function <f> defined in M a x K such that ||0||2,p,^,<7 < 1- Then, 

|2VCiV0 + 0ACi| < C e -^(l+r( ra ))-l^|| 2w . 
We also have that ||S(w)|| PiM-(T < Co? , hence 

|(l-C)S(w)| <Ce-^(l + r (ax))-" 

and 

||<?i|U < Ce"< 

Let consider the set 

A = {^eX / \\i,\\ x <Ae^i}, 

for a large number A > 0. Since 

I K{ih) - K(i> 2 ) I < C(l - Ci) sup |iV>! + (1 - t)^, + (201 - 02| , 

te(o,i) 



we find that 



||A'(^ 1 )-A'(^ 2 )|| 00)|1 < Ce-^||Vi-02|| 



while ||i^(0)||oo,/i < Ce"* 7 " . It follows that the right hand side of equation (|7.5|) 
defines a contraction mapping of A, and hence a unique solution ip — ^(4>) G A 
exists, provided that the number A in the definition of A is taken sufficiently large 
and ||0||2,p,/i,<r < 1- In addition, it is direct to check the Lipschitz dependence of ^ 
(EH) on ||Vlkp,^ < 1. □ 

Thus, we replace replace ip = *&(4>) into the equation (|4.14|) of the gluing system 
(|4.12p - (|4.14p and get the (nonlocal) problem, 

dud) + A y , M > = -S(ui) - M(0) in M a x R (7.6) 

where 

M(0) := B(0) + - j»]0 + Cl(/'(«l) - + Cl#TO) + 0), 

V v ' " v ' V , ' 

NlO) N 2 (» N 3 (0) 

(7.7) 

which is what we concentrate in solving next. 
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7.2. Proof of Proposition l4~2l We recall from M.2\ that Proposition 14.21 refers 
to solving the projected problem 

d tt (f> + A vMa (f> = -§(ux) - N(0) + c(y)w'(t) in M a x R, 

f (7-8) 
/ (f>(y, t) w'{t) dt = 0, for all y e M a , 

and then adjust h\ so that c(y) = 0. Let <j) — T(g) be the linear operator providing 
the solution in Proposition 14.11 Then Problem (|7.8j) can be reformulated as the 
fixed point problem 

<t> = Tf-SM-NMJ-Tfo), ||0|| 2 „<1 (7.9) 
which is equivalent to 

<t> = T(-S( Ul )+a 2 Ahiw' -U(ct>)), Hh m <l, (7.10) 

since the term added has the form p{y)w' which thus adds up to c(y)w'. The reason 
to absorb this term is that because of assumption (14. 7|) . ||a 2 A/ii iy'|| P) 4 )CT = 0(a 3 ~p) 
while the remainder has a priori size slightly smaller, 0(a 3 ). 

7.2.1. Lipschitz character o/N. We will solve Problem (|7.10[) using contraction map- 
ping principle, so that we need to give account of a suitable Lipschitz property for 
the operator T. We claim the following. 

Claim. We have that for a certain constant C > possibly depending on K, in 
( |^. 7[ ) but independent of a > 0, such that for any <fii, i>i with 

||N(0i) - N(0 2 )|U M +i,a < Ca\\(fn- hh, P ,^ (7.11) 
where the operator N is defined in | | 7.7| . 

We study the Lipschitz character of the operator N through analyzing each of 
its components. Let us start with N±. This is a second order linear operator with 
coefficients of order a plus a decay of order at least (^(r" 1 ). We recall that B = Q2.B 
where in coordinates 

B = (f '(ui) - /' ») - a 2 [(t + hi)\A\ 2 + & M hi]dt -2aa%djhd u + 

a(t + h) [aljdij - a al^djhdu + dihd jt ) + a(b\di - ab\d l hd t ) ) ] + 

a 3 {t + hfb\dt + a 2 [a% + u{t + hy^dihdjhdtt (7.12) 
where, we recall, 

^ = 0(r- 2 ), a^ = 0(r~ 2 ), b\ = 0(r- 3 ), bf = 0(r- 6 ), 

/'(ui) -/'(«;) =0(a 2 r- a e-'l*l) 8^ = 0^), \A\ 2 = 0(r^) . 
We claim that 

IINiMIUh-i,* < C«ll^ll2 >P> ^. (7.13) 
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The only term of iVi(0) that requires a bit more attention is a 2 (Ahi)(ay)dt<j) ■ We 
have 

f \a 2 (Ah 1 )(az)d t (j)\ p dV a (z)dT < 

JB((y,t),l) 

Ca 2 P||^|| L oo (B(( ^ ) . 1) (l+r a (y))- 4 f+ 4 f \(l+r a (z) ) 4 ~f (Ah^az)^ \dV a (z) < 

JB((y,t),l) 

and hence in particular for p > 3, 

\\a (A/ii)(ay)e^||p )(tl+ 2 )0 - < C a * ||^i||* Iwkp,/^ < C a v\\(t>\\2,p,w 
Let us consider now functions 0; with 

||<&||2,p,/i,<T < 1, i = 1,2. 

Now, according to Lemma [47T1 we get that 

||N 2 (0i) - N 2 («/» 2 )|| P)ft(r < Ce-'^Wfa - <h\\vw ■ (7-14) 
Finally, we also have that 

1^(00 -N 3 (0 2 )| < 
CCl Sup |t(*(0l) + 0l) + (l-t)(*(02) + 02)| -02| + |*(^)-*(02)|], 

te(o,i) 

hence 

||N 3 (0i) - N 3 (0 2 )||p,2 M , ff < C ( ||0l||oo. M ,<T + ||02||oo. M ,<t +e _<T «) ||01 - 2 || OO . AI . CT . 

(7.15) 

From (|7.13[) . (|7.14[) and (|7.15[) . inequality (|7.1ip follows. The proof of the claim is 
concluded. 



7.2.2. Conclusion of the proof of Proposition \4-2\ The first observation is that 
choosing /i < 3, we get 

\\S( Ul ) + a 2 Ah 1 w'\\ p ^ a < Ca 3 . (7.16) 

Let us assume now that 0i , 02 G B a where 

B a = {0 / Uh,p,„,c < Ka 3 } 

where if is a constant to be chosen. Then we observe that for small a 

||N(0)|| p , M +i, CT < Ca 4 , for all G B a , 

where C is independent of K. Then, from relations 1|7.16[) - (|7.15[) we see that if K 
is fixed large enough independent of a, then the right hand side of equation (|7.5[) 
defines an operator that applies B a into itself, which is also a contraction mapping 
of B a endowed with the norm || Hp.^o-, provided that /i < 3. We conclude, from 
contraction mapping principle, the existence of as required. 

The Lipschitz dependence (I4.27P is a consequence of series of lengthy but straight- 
forward considerations of the Lipschitz character in hi of the operator in the right 
hand side of equation (|7.5[) for the norm || ||* defined in (|4.34[) . Let us recall expres- 
sion (|7.12[) for the operator B, and consider as an example, two terms that depend 
linearly on hf. 

A(hi,4>) := aa°j djhid U (f> ■ 

Then 

\A(h u <t>)\ < Ca\djhi\\d it <j>. 



MANUEL DEL PINO, MICHAL KOWALCZYK, AND JUNCHENG WEI 



Hence 

\\A(hi,<j>)\\ PlM .2, a < Ca||(l + r2)^/ii|U||a it 0lU, CT < Ca 4 ||/n||* Uh, P ^- 
Similarly, for A((f>, hi) = q; 2 Aa//ii dt<j> we have 

\A{ct>,h 1 )\ < Ca 2 |A M /i 1 (ay)|(l+r Q )-' 1 e- CT l t l||^|| 2w . 

Hence 

\\a 2 A M hi 9*01^,^+2,0- < Ca^J \\hi\\* \\<f>\\2,p,n,a- 

We should take into account that some terms involve nonlinear, however mild de- 
pendence, in hi. We recall for instance that a\j = a}j(ay, a(t+ho + hi)). Examining 
the rest of the terms involved we find that the whole operator N produces a depen- 
dence on hi which is Lipschitz with small constant, and gaining decay in r a , 

||N(/ii,0) - N(/i 2 , <f>)}\ P ,„+i,a < Ca 2 \\hi - h 2 \\* Uh, P ,w (7.17) 

Now, in the error term 

TZ = -5(iti) + a 2 A/iiw', 

we have that 

\\K{hi) - TZ{h 2 )\\ p , 3 ,a <Ca 2 \\hi- h 2 \U • (7.18) 
To see this, again we go term by term in expansion (|4. 15|) . For instance the linear 
term a 2 afjdihodjhi w" . We have 

\a 2 a^OMdM < C a 2 (I + r a )- 3 e-°W \\hi\U 

so that 

\\a 2 a°- dihodjtuWp^a < C a 2 \\hi\\*. 

the remaining terms are checked similarly. 

Combining estimates (|7.17|) , (|7. 18[) and the fixed point characterization (|7.5|) we 
obtain the desired Lipschitz dependence (14.27| of $. This concludes the proof. □ 

8. The reduced problem: proof of Proposition ^. J\ 

In this section we prove Proposition 14.41 based on the linear theory provided by 
Proposition ^. 31 Thus, we want to solve the problem 

J 

J(hi)=A M h 1 +hi\A\ 2 =G(hi)+J2— ^-z, inM, (8.1) 

i—1 

[ dV = for all i — 1 ■ ■ ■ J 

where the linearly independent Jacobi fields z% will be chosen in (|9.1| and (|9.2jl of 
§8, and G = G1+G2 was defined in (|4.29ll . (|4.30| . We will use contraction mapping 
principle to determine the existence of a unique solution hi for which constraint 
(|4.7|) . namely 

||^i||*:=||/ii||L-(M) + ||(l + r 2 )^i||L~(M) + ||^ 2 ^i|| P ,4-| < (8-2) 

is satisfied after fixing K, sufficiently large. 

We need to analyze the size of the operator G, for which the crucial step is the 
following estimate. 



THE ALLEN CAHN EQUATION AND MINIMAL SURFACES IN R 3 



39 



Lemma 8.1. Let tp(y,t) be a function defined in M a x R such that 

II^Hjw, := sup e a ^(l + r^)\\iP\\ LPiB{iVtthl) < +00 

(y,t)eM a xR 

for a, /i > 0. The function defined in M as 

q(y) := [ iP(y/a,t)w'{t)dt 
Jr 

satisfies 



provided that 

In particular, for any t > 0, 
and 



\q\\ P ,a < c 



u > — ha. 

p 



IMI P ,2-f-r < C\m P ,2,c 



(8.3) 



||g|| Pl4 -i < C||V|| P ,4, CT . 

Proof. We have that for \y\ > R 



f \yr [ i>{y/a,t)w' 

J\y\>Ra JWL 



Now 



(t)dt 



dV < C 



[ w'(t)dt [ 
Jm J\v\>Rq 



(8.4) 
(8.5) 

\yr k%AM)PW. 



and 



f \ y r\i>{y/a,t)\vdV = a a ?+ 2 f \y\" 

J\y\>Ro J\y\>Ro/a 

[ \vrm,t)\>w a < c Yl * ap [ 

J\y\>Ro/a i>[Ro/a] J '*<lwl<« 



\i>(y,t)\ p dV a 



\Hy,t)\ p dv a 



Now, i < \y\ < i + 1 is contained in 0(i) balls with radius one centered at points of 
the annulus, hence 

Ji<\y\<i+1 

< Ce-^MW^ I (ltr„)-'W a 

Ji<\y\<i+\ 



Ji<\y\<i+1 



Then we find 



|a g l|P 



LP(\y\>R ) - ° a imi P)M * 

»>[flo/a] 



The sum converges if /1 > ^ + a and in this case 



so that 



|y| a 9lU P (| y |>K„) < C\\i/>\\ Plll . 
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Now, for the inner part \y\ < Rq in M, the weights play no role. We have 
\i/)(j/a,t)\*dV = a 2 f \^(y,t)\ p dV a < 

v\<Ro J\y\<R /a 

Co? ! mV,t)\ P dV a <Ca 2 W\\^e-°*>W £ i 

i<R / a Jl <y\ <l+1 i<R /a 

Hence if fi > | + a we finally get 



and the proof of (|8.3[) is concluded. Letting (/i, a) = (2, 2— | — r), (/i, a) = (4, 4 — ~) 



lkllp,a < CIIV'IIp.m 

4 - 
P' 

respectively in ((53)) . we obtain (JHHJ) and f53]) . □ 

Let us apply this result to ip(y, t) = N($(/ii) ) to estimate the size of the operator 
G 2 in (jlUO) . For = $(/n) we have that 

G 2 {hx){y) -.^cZ'a- 2 [ N(4>)(y/a,t)w'dt 

JR 

satisfies 

||G 2 (/n)|| Pi4 _| <Ca- 2 ||N(0)|| p , 4iff < C« 2 . 
On the other hand, we have that, similarly, for (f>i — &(hi), I = 1,2, 

\\G 2 (h 1 )-G 2 (h 2 )\\ 4 _, < Ca- 2 ||N(^,^)-N(0 2 , h 2 )\\ pA . a . 

' p 

Now, 

||H(0i,^i) - N(0x, /i 2 )||p, 4 , - < Ca 2 ||/ii - M*||</>i|kp,3,<T, < Ca 5 ||/n - /i 2 ||*, 
according to inequality (|7.17[) , and 

Picpuhi) - N(0 2 , fti)|| P) 4, CT < Ca 2 ||0! - 2 ||p,3,a < Ca 4 ||/H - /i 2 ||* . 
We conclude then that 

||G 2 (/ii)-G 2 (/ l2 )|L 4 _4 < Ca 2 \\h 1 -h 2 \U ■ 
In addition, we also have that 

||G 2 (0)|| p . 4 _4 <Ca 2 . 

for some C > possibly dependent of JC. On the other hand, it is similarly checked 
that the remaining small operator Gi(/ii) in (14.29[) satisfies 

||Gi(/n)-Gi(/i 2 )|L 4 _4 < Cia\\hi-ha\\* . 
A simple but crucial observation we make is that 

c*Gi(0) = adihodjho / Ci(* + M a ij W c?t + a~ 2 / C4 Ri(y, t, 0, ) w' dt 

JR JR 

so that for a constant C 2 independent of /C in (|8.2|) we have 

||Gi(0)|| p , 4 _4 <G 2 a. 

In all we have that the operator G{h\) has an 0(a) Lipschitz constant, and in 
addition satisfies 

||G(0)|| Pl 4-i <2G 2 a. 
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Let h = T(g) be the linear operator defined by Proposition 14.31 Then we consider 
the problem (|8.ip written as the fixed point problem 

hi = T(G(hi)), \\h\\* < ICa. (8.6) 

We have 

\\nG(hi))\\. <||T||||G(0)|| Pi4 _ |+ Ca||^||,. 

Hence fixing JC > 2C*2 ||T r || , we find that for all a sufficiently small, the operator 
TG is a contraction mapping of the ball ||/i||* < ICa into itself. We thus have the 
existence of a unique solution of the fixed problem (|8.6p . namely a unique solution hi 
to problem (|8.1|) satisfying (|8.2p and the proof of Proposition 14.41 is concluded. □ 

9. Conclusion of the proof of Theorem]^ 
We denote in what follows 



r{x) = Jx\ + f = -(xi,x 2 ,0), 9 =-(-x 2 ,Xi,0) . 

V J* 



We consider the four Jacobi fields associated to rigid motions, Z\, . . . , Z4 intro- 
duced in (|1.14p . Let J be the number of bounded, linearly independent Jacobi 
fields of J . By our assumption and the asymptotic expansion of the ends (|1.12[) . 
3 < J < 4. (Note that when M is a catenoid, 2:4 = and J = 3.) Let us choose 

4 

'= z2 d ji z oi,3 : = 1,-, J (9-1) 
{=1 

be normalized such that 

q(y)ziZj = 0, for i ± j, q(y)zf = 1, i, j = I, ■ ■ ■ , J. (9.2) 

M JM 

In what follows we fix the function q as 

i(y) ■■= 1x n 4 ■ (9.3) 

So far we have built, for certain constants c 2 ; a solution u of equation (|4. 36[) . 
namely 

J 

Au + f(u) = ^2 CiZi(ay)w' (t)q(ay)( 2 

where u, defined in (|4.35|) satisfies the following properties 

u(x)=w{t)+<f>(y,t) (9.4) 
near the manifold, meaning this x — y + (t + h(ay) ) ^(cny) with 

y€Af Q , |t| < -+ 7 log(2 + r(ay)). 

a 

The function </> satisfies in this region the estimate 

|0| + |V0|<g a 2 1 e-^l. (9.5) 
1 + r z {ay) 

Moreover, we have the validity of the global estimate 

C 



|Vu(x)| < 1 
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We introduce the functions 

Zi(x) — d Xi u(x), i — 1,2, 3, Zi(x) — —ax2d X2 u + ax\d X2 u . 
From the expansion (|9 .4(1 we see that 

Vu(x) = w'(t)Vt + V0. 

Now, t = z — h(ay) where z designates normal coordinate to M a . Since Vz = v 
v(ay) we then get 

Vt = v{ay) — aVh(ay). 
Let us recall that h satisfies h = (— l) k (3k logr + 0(1) along the fc-th end, and 

V/i= (-l) fc ^f + C>(r- 2 ). 
r 

From estimate (19.51) we we find that 



Vu{x) = w'itMv - a{-l) k —f) + 0{ar' 2 e~ a ^). (9.6) 
From here we get that near the manifold, 

Zi{x) =w'(t)(z l (ay)-a(-l) k ^re l ) + 0(ar- 2 e-^), i = 1,2,3, (9.7) 



r a 

Z A {x) = w'{t) zoi(ay) + Oiar^ 1 e^ 1 * 1 ). (9.8) 

Using the characterization (|4.36p of the solution u and barriers (in exactly the 
same way as in Lemma \1 1 .41 below which estimates eigenfunctions of the linearized 
operator), we find the following estimate for r a {x) > Rq; 

m 

fc=i 

We claim that 

(Am + f{u))Z t {x)dx = for all i = l,...,4 (9.10) 
so that 

Cj / q(ax)zj(ay)w'(t)Zi(x)(2dx = for all i = l,...,4. (9.11) 

Let us accept this fact for the moment. Let us observe that from estimates (19.7[) 
and EH), 



a 2 / q(ax)zj(ay)w' (t) duZi(x) (2 = / w'{t) 2 dt I qZjZidV + o(l) 

with o(l) is small with a. Since the functions z, are linearly independent on any 
open set because they solve an homogeneous elliptic PDE, we conclude that the 
matrix with the above coefficients is invcrtiblc. Hence from (|9.1ip and (|9.2[) . all 
Cj's are necessarily zero. We have thus found a solution to the Allen Cahn equation 
(II. ip with the properties required in Theorem [T] 
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It remains to prove identities (|9.10| . The idea is to use the invariance of A + f(u) 
under rigid translations and rotations. This type of Pohozaev identity argument 
has been used in a number of places, see for instance [15] . 

In order to prove that the identity (]9 . 10[) holds for i — 3, we consider a large 
number R >> ^ and the infinite cylinder 

C R = {x I x\ + xj < R 2 }. 

Since in Cr the quantities involved in the integration approach zero at exponential 
rate as \xa\ — > +oo uniformly in (xi, X2), we have that 

(Au + f(u))d X3 u- [ Vu-rd X3 u= [ a x3 (F( U )-i|Vu| 2 ) = 0. 

C R JdCn JCr 1 



We claim that 



-R— >+oo 



lim / Vtt • f d X3 u = 0. 



dC R 



Using estimate (|9.6[) we have that near the manifold, 

d X3 uVu(x) ■ f = w'(t) 2 ((u - a{~l) k —f) • rV 3 + Ofae - "'*! \). 

r a r z 

Let us consider the k-ih end, which for large r is expanded as 

%3 = F k:a (x 1 ,x 2 ) = a~ 1 (a k \ogar + b k + O^ 1 )) 

so that 

{-Ifv = — f====== ( VFfc Q , -1) = ^- - e 3 + 0(r- 2 ) . (9.12) 

|V^., Q | 2 a r 

Then on the portion of Cr near this end we have that 

[y — a( — l) — r) ■ ris 3 = -a h OiR ). (9.13) 

r a R 

In addition, also, for x 2 + x\ = R 2 we have the expansion 

t= (x 3 - F k<a {x 1 ,x 2 ) - falogar + 0(1))(1 + 0(R~ 2 )) 

with the same order valid after differentiation in X3, uniformly in such (0:1,2:2). 
Let us choose p — ^\ogR for a large, fixed 7. Observe that on 8Cr the distance 
between ends is greater than 2p whenever a is sufficiently small. We get, 

rF k , a (xi,X2)+/3k logctr+p rco 

/ w / (t) 2 dx 3 = w'{t) 2 dt + 0{Rr 2 ). 

J -Ffe, « (ail, ^2)+/3 fc log ai — p J ~oc 

Because of estimate (19.9[) we conclude, fixing appropriately 7, that 



As a conclusion 



d X3 uVu{x) ■fdx i = 0(R~ 2 ) . 
f} k {\x3-Fk, a \>p} 



d X3 uVu-rdx 3 = ^^2{a k + a(3 k ) / w' {tf dt + 0{R- 2 ) 

and hence 

/ 3 l3 uV«(j)-r= VK + aftl + Or 1 ). 
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But Y^k=i °fe — Ylk=i Pk — and hence (|9.10|) for i = 3 follows after letting 
R — ► oo. 

Let us prove the identity for i = 2. We need to carry out now the integration 
against d X2 u. In this case we get 

/ (Au + f(u))d X2 u= I Vu-r8 X2 u+ [ d X2 (F{u)~l\Vu\ 2 ). 

JCr JOCr JCr 1 

We have that 

d X2 (F(u)-l\Vu\ 2 ) = ( (F(u)-l\\7u\ 2 )n 2 



Cr 2 JdC R 2 

where n 2 = x 2 /r. Now, near the ends estimate (|9.6[) yields 

\Vu\ 2 = \w'(t)\ 2 +0(e-"W±) 
and arguing as before, we get 

/oo 
\w'(t)\ 2 dt + 0(R- 2 ). 
-oo 

Hence 

/ \\7u\ 2 n 2 =m [ \w'{t)\ 2 dt [ n 2 + 0{R- 1 ). 

JOCr J-oo J[r=R] 

Since Jr r — R i n 2 = we conclude that 

lim / \Wu\ 2 n 2 =0. 

R^+oc J gCli 

In a similar way we get 

lim / F(u) n 2 = 0. 

R^+oo Jq Cr 

Since near the ends we have 

d X2 u = w'(t)(v 2 - a (-l) k ^re 2 ) + 0(ar- 2 e -^) 

and from (|9.12ll v 2 = 0(R^ 1 ), completing the computation as previously done 
yields 

/ Vu -fd X2 u = 0{R- 1 ). 

JdC R 

As a conclusion of the previous estimates, letting R —> +oo we finally find the 
validity of (|9.10|) for i — 2. Of course the same argument holds for i = 1. 

Finally, for i = 4 it is convenient to compute the integral over Cr using cylin- 
drical coordinates. Let us write u = u(r, 0, z). Then 

(Au + f(u)) (x 2 d Xl u - xid Xl u) = 

>Cr 

r2ir pR poo 

[uzz + r~ 1 (ru r ) r + f(u)]ugr d8 dr dz = 

>0 JO J -oo 

f27r pR. poo poo p2ir 



'0 JO J -oo 



/OO pOO pZTT 

de[u 2 z +u 2 r ~2F(u)]rdddrdz+R / u r u e (R,9, z) d9 dz 
-OO J —OO JO 
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+ / U r Ue ■ 
JdC R 

On the other hand, on the portion of BCr near the ends we have 

u r u e = w'{t) 2 R{v ■ f){v ■ 9) + 0(R- 2 e~ a W). 
From (f97T2|) we find 

(v ■ r)(v ■ 0) = 0{R- 3 ), 

hence 

u r ug = w'(t) 2 0(R- 2 ) + 0{R- 2 e- aW ) 

and finally 

u r u e = OiR- 1 ). 



ldC R 

Letting R — » +00 we obtain relation (|9 . 1 0[) for i = 4. The proof is concluded. □ 



10. Negative eigenvalues and their eigenf unctions for the Jacobi operator 

For the proof of Theorem [5] we need to translate the information on the index of 
the minimal surface M into spectral features of the Jacobi operator. Since M has 
finite total curvature, the index i(M) of the minimal surface M is finite. We will 
translate this information into an eigenvalue problem for the operator J . Let 

Q(k,k) := f \Vk\ 2 dV- [ \A\ 2 k 2 dV. 
Jm Jm 

The number i(M) is, by definition, the largest dimension for a vector space E of 
compactly supported smooth functions in M such that 

Q(z,z)<0 for all zeE\{0}. 

The number i(M) when finite has the following convenient characterization, 
whose proof is straightforward. In what follows we fix the function q as 

p(y) ■= tt^i • do-14) 

Let us consider for a large number R, the region 

M R = {yeM I r(y) < R} 

and the eigenvalue problem 

A M k+\A\ 2 k + \p{y)k = mAI fl , (10.15) 
k = on dM R . 

Let 771r(p) denote the number of negative eigenvalues (counting multiplicities) for 
this problem. Then we have 

i(M) = sup m R (p) . (10.16) 

R>0 

Let us also consider the eigenvalue problem in entire space 

A M k + \A\ 2 k + Xp(y)k = in M, k^L co {M). (10.17) 
We will prove the following result. 

Lemma 10.1. Problem ( | 10. 17^ has exactly i(M) negative eigenvalues, counting 
multiplicities. 
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10.0.3. A priori estimates in M . For the proof of Lemma 110.11 and for later 
purposes, it is useful to have a priori estimates uniform in large R > for the linear 
problem 

A M k+\A\ 2 k- 1P {y)k = f mM R , (10.18) 
k = on dM R . 

We have the following result. 

Lemma 10.2. Let p > 1, a > 0. Then for Rq > large enough and fixed and 
7o > ; there exist a C > such that for all R > Rq + 1, < 7 < 70, any f , and 
any solution k of problem (j 10. 18$ . we have that 
(«) //I/L4-A < +00 then 

||fc||oc < C[||/llp,4-f + l|fc|U-(|»|<3flo) ] ■ (10-19) 

(b) #ll/L,2- = -<7 < +°°; then > 

\\D 2 k\\ p ^__ a + || Dk \\ p>l _z_ a < C [ ||/|| P)2 _ 2 _ CT + ||fc||oo ] ■ (10.20) 

f 5 p f 1 p rip 

If p > 2, we have in addition 

IKi + M) 1 -^!!^ <c[||/|| Pi2 _2_ ff + WkU]. (10.21) 

Vi p 

Proof. Let us consider the equation in M 

A M Tp + \A\ 2 iP = -\f\x\y\ <R , \y\>Ro, (10.22) 

ip(y) = 0, \y\ = Ro. (10.23) 
For a large and fixed Rq , solving this problem amounts to doing it on each separate 
end. As in Lemma l6.1[ after a Kelvin's transform the problem reduces in each end 
to solving in a ball in R 2 an equation of the form 

AV> + 0(\y\ 2 )D 2 4, + 0(\y\)D{j + 0(1$ = -|/|X| y |>A, Ixl < ^ , 

Rq 

$(y)=o, |y| = i • 

Enlarging i? if necessary, this problem has a unique solution, which is also positive. 
This produces a bounded, positive solution ip of (|10.22[) - (|10.23p with 

Moo < c\\f\\ pA _, . 

On the other hand, on this end the Jacobi field Z3 = v ■ e% can be taken positive 
with Z3 > 1 on |y| > Rq. Thus the function ip + ||fc||z, oo (|j/|=Ro)^3 is a positive, 
bounded supersolution for the problem (|10.18[) in this end, where |y| > Rq, and 
estimate ()10. 19|) then readily follows. 

Let us prove now estimate (|10.20[) . Fix a large number Rq > and another 
number R >> Rq. Consider also a large p > with ip < R. On a given end we 
parameterize with Euclidean coordinates y £ K 2 and get that the equation satisfied 
by k = fc(y) reads 

Ak + 0(\y\- 2 )D 2 k + 0(\y\- 3 )Dk + 0(|y|- 4 )fc = /, Rq < |y| < R. 

Consider the function k p (z) — k{pz) wherever it is defined. Then 

Ak p + 0(p- 2 \z\- 2 )D 2 k p + 0(p- 2 \z\- 3 )Dk p + 0{p- 2 \z\- A )k p = f p 
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where f p (z) = p 2 f(pz). Then interior elliptic estimates (see Theorem 9.11 of [14] ) 
yield the existence of a constant C = C(p) such that for any sufficiently large p 

\\^p\\LP{l<\z\<2) + l|-D 2 fcpllLP(l<|2|<2) < C ( l|fcplli«>(£<|*|<3) + ll/plliP(i<|z|<3) )• 

(10.24) 

Now, 

II/pIIl,(K N <3) =P 2P [ \f(pz)\ p dz< 
C P pa [ \ P z\ 2p ' 2 -^\f( P z)\ p p 2 dz = C(?° [ \y\ 2p ~ 2 - pa \f(j)\ p dy. 

J(i<\z\<3) J(%<\y\<3p) 

Similarly 

IP 2 MW|*i<2) > cr ! \y\ 2p - 2 - pa \D 2 k{ y )\vd Y . 

J( P <\y\<2p) 

Thus 

f lyl 2 P -2- P a lD 2 k{yrdY < c f |y| 2p - 2 - pCT |/(y)l p rfy+p^l|fc|IL 

J(p<\y\<2p) J(£<| y |< 4 p) 

Take p = pj = IK Then 



f \y\ 2p ~ 2 ~ P °\D 2 k(y)\rdy < 



'(pj<|y|<Pj+i) 

C f \y\ 2p ~ 2 - pa \f{y)\ p dy +2-^||fc|| p oo- 

Then, adding up these relations wherever they are defined, taking in addition into 
account boundary elliptic estimates which give that for p = -j, 

l|-D 2 MU p (l<M<2) < C ( ||fcp||L° c (i<| 2 |<2) + ll/pllLP(i<|z|<2) ) J 

plus a local elliptic estimate in a bounded region, we obtain that for some C > 
independent of R, 

ll# 2 fc||p,2-f -a < CdlfcHoc + ll/IUa-l-J. 

The corresponding estimate for the gradient follows immediately from (|10.24p . 
We have proven (|10.20|) . If p > 2 we can use Sobolev's embedding to include 
||-DfcpliL°=(i<|z|<2) on the left hand side of (j!0.24|) . and estimate (|10.21|) follows. 
The proof is complete. □ 



10.0.4. Proof of Lemma \10.1{ We will prove first that problem (|10.17|) has at 
least i(M) linearly independent eigenfunctions associated to negative eigenvalues 
in L°°(M). For all R > sufficiently large, problem (|10. 15[) has n = i(M) linearly 
independent eigenfunctions ki,R, ■ ■ . , k ni R associated to negative eigenvalues 

< ^2M < ' ' ' < KiM < . 

Through the min-max characterization of these eigenvalues, we see that they can 
be chosen to define decreasing functions of R. On the other hand, Ai^ must be 
bounded below. Indeed, for a sufficiently large 7 > we have that 

\A\ 2 - 7p < in M 



48 



MANUEL DEL PINO, MICHAL KOWALCZYK, AND JUNCHENG WEI 



and by maximum principle we must have Xi,r > —7. The eigenfunctions can be 
chosen orthogonal in the sense that 

/ pk lR k jR dV = for all i^j. (10.25) 
Jm r 

Let us assume that 1 1 fc» , j£ 1 1 00 = 1- Then the a priori estimate in Lemma 110.21 imply 
that, passing to a subsequence in R — > +00, we may assume that 

Xi,R I Aj < 0, k i<R (y) -> ki(y), 

uniformly on compact subsets of M, where ki ^ is a bounded eigenfunction of 
(|10.17p associated to the negative eigenvalue A^. Moreover, relations (|10.25[) pass 
to the limit and yield 

/ phkjdV = for all i^j. (10.26) 

JM 

Thus, problem (|10. 17|) has at least n — i(M) negative eigenvalues. Let us assume 
there is a further bounded eigenfunction linearly independent of ki, . . . , k n , 

say with 

/ pk t k n+1 dV = for all i = l,...,n, (10.27) 

JM 

associated to a negative eigenvalue A„+i. Then the a priori estimate of Lemma l6.ll 
implies that 

||(1 + r 2 )Vfc„+i|| <+oo. 

The same of course holds for the remaining fej's. It follows that 

Q(fc,fc)<0 for all k G span{fci, . . . , k n +i} \ {0}. 

However, again since Vfcj decays fast, the same relation above will hold true for 
the fci's replaced by suitable smooth truncations far away from the origin. This 
implies, by definition, i(M) > n + 1 and we have have reached a contradiction. The 
proof is concluded. □ 



11. The proof of Theorem^ 

In this section we will prove that the Morse index m(u a ) of the solution we have 
built in Theorem [1] coincides with the index of the surface M, as stated in Theorem 
[5] We recall that this number is defined as the supremum of all dimensions of vector 
spaces E of compactly supported smooth functions for which 

Q(i/j,i/f) = / |V^| 2 -/'(u Q )V 2 <0 for all £ E \ {0}. 



We provide next a more convenient characterization of this number, analogous to 
that for the Jacobi operator of ^ 101 Let us consider a smooth function p(x) defined 
in R 3 such that 

p(ax) = — — prj if x = y + (t + h{ay))v(ay) eAfs, 

J- + r a (y) 

and such that for positive numbers a, b, 

< p(ax) < p — rj for all x = ( x' , x 3 ) e K 3 . 



1 + lax'l 4 ~ ^ v ' ~ 1+ lax' 
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For each R > 0, we consider the eigenvalue problem in the cylinder 
C R = {(x',x 3 ) I \x'\ < Rot' 1 , \x 3 \ < Ra- 1 )}, 

A<t> + f'(u a )<l> + \p(ax)(t> = inC R , (11.1) 

= on 8Cr . 

We also consider the problem in entire space 

A0 + /'(u Q )0 + Xp(ax)cj) = inM 3 , 0eL°°(IR 3 ). (11.2) 

Let mR{u a ) be the number of negative eigenvalues A (counting multiplicities) of 
this Problem (jll.ip . Then we readily check that 

m{u a ) = sup m R (u a ). 

R>0 

On the other hand, we have seen in flOlthat the index i(M) of the minimal surface 
can be characterized as the number of linearly independent eigenfunctions associ- 
ated to negative eigenvalues of the problem 

Az + \A\ 2 z + Xp(y)z = in M, zeL°°(M), (11.3) 

which corresponds to the maximal dimension of the negative subspace in L°°(M) 
for the quadratic form 

Q{z,z) = [ \V M z\ 2 -\A\ 2 z 2 dV. 

J M 

We shall prove in this section that m(u a ) = i(M) for any sufficiently small a. 

The idea of the proof is to put in correspondence eigenfunctions for negative 
eigenvalues of problem (jll.ip for large R with those of problem (| 11 . 3(1 . This cor- 
respondence comes roughly as follows. If z is such an eigenfunction for problem 
(|11.3p then the function defined near M a as 

k(y)w'(t), k(y)=z(ay) (11.4) 

defines after truncation a negative direction for the quadratic form Q on any large 
ball. Reciprocally, an eigenfunction for negative eigenvalue of problem (jll.ip will 
look for any sufficiently small a and all large R like a function of the form (jll.4p . 
In the following two lemmas we clarify the action of the operator L on functions of 
this type, and the corresponding connection at the level of the quadratic forms Q 
and Q. 

Lemma 11.1. Let k{y) be a function of class C 2 defined in some open subset V of 
M a . Let us consider the function v(x) defined for x £ Ms, y S V as 

v(x) = v(y,t) := k(y)w'(t), y e V, \t + h 1 {ay)\ < p a {y) 

where p a is the function in the definition ofJ\f a , (\3. 11$ . Then L(v) := A x v+f'(u a )v 
can be expanded as in (| 11.8$ below. Besides we have 

[ L{v)w'dt = {A Ma k + a 2 \A\ 2 k + ahalfdtjk) [ w' 2 dt 

J\t+hl\<p a JR 

+ 0{a 2 r- 2 ) dijk + 0(a 2 r- 3 ) 8,k + 0(a 3 r - 4 ) k . (11.5) 

Here 

a lf = a lf(ay) = 0{r~ 2 ) . 
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The same conclusions hold for the function 

v(x) = v(y,t) := k(y)w'(t)r)s{y,t), y G V, \t + h 1 {ay)\ < p a (y) 
where the cut-off function r)$ is defined in (jff. 12\ . 
Proof. Let us recall that 

A x = d tt + A Ma - a 2 [(t + h)\A\ 2 + A M h}d t - 2a o° djhdtt + 
a(t + h) [aljdij — 2a ajj djhdu + a{h\di — abjdihdt) ) ] + 
a 3 (t + h) 2 b\d t + a 2 [a% + a(t + h)al j )]d i hd j hd tt . 
Hence, using Lemma 12.21 in the appendix we get 

A x v + f'(u a )v = k(w"' + f'(w)w') + [/'(«„) - f'(w)} kw' + 
w' A Ma k - a 2 [(t + h 1 )\A\ 2 + A M h!]kw" - 2aa° ij d j hd i kw" + 



(t + h) [aljdijkw' — a a\Adjhdik + dihdjk) w" + a(b\dik w' — ab\dihw") ] + 



n 



{t + h) 2 b\kw" + a 2 [a°- + a(t + h)alA}d l hd j hkw"' . (11.6) 



We can expand 

alj = a]j(ay, 0) + a{t + h) af^ay, a(t + h)) =: a\f + a(t + h)a 2 j, 
with a 2 j = 0(r~ 2 ), and similarly 

b) = b]{ay, 0) + a{t + h) b 2 (ay, a{t + h)) =: b 1 / + a(t + h)b 2 , 
with b 2 = 0(r~ 3 ). On the other hand, let us recall that 

u a - w = fa + 0(a 3 r- 4 e- CT l*l) 
where fa is given by (|3.2[) . 

fa(y,t) = a 2 \A(ay)\ 2 M-t) ~ oi 2 a%d i h d j h {ay)fa{t) (11.7) 
and V'Oj fa decay exponentially as |t| — * +oo. Hence 

[f'(u a )-f'(w)]w' = f"{w)faw' + 0(a 3 e-^r- 4 ). 
Using these considerations and expression (|11.6|) . we can write, 

Q := A x v + f'{u a )v = 
A Ma kw' - a 2 \A\ 2 ktw" + a 2 a^dihodjho kw'" + aha 1 l fd tj kw' + f"{w) fa kw 1 



-w 



Qi 

aa^Adjhdik + dihdjk) + a 2 kAnjh\ + a 2 ha\f (djhdik + dihdjk) 



atw 



a ij + a V 
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+ a 2 (t + hfalfiijkw' + a 2 (t + h)a 2 Ad 3 hdik + d l hd j k) w" + 0(a 



(11.8) 

The precise meaning of the remainder Qq is 

Q 6 = 0(a 3 e-^r~ 2 )d t3 k + 0(a 3 e~^r~ 3 )d 3 k. 
We will integrate the above relation against w'(t) in the region \t + hi(ay)\ < p a (y). 

Let us observe that the terms Qi for i — 1, . . . , 4 are in reality defined for all t 
and that 

f Qiiv'dt = [ Q lW ' dt +0(a 3 r- i ) (11.9) 

J \t+h-i_\<p a JM 

where the remainder means 

0(a 3 r~ 4 ) := 0{a 3 r- 4 ) d tJ k + 0(a 3 r~ A ) d t k + 0(a 3 r - 4 ) k . 
Let us observe that 

(Qs + Qi)w'dt = . (11.10) 



On the other hand, since 

J tw"w'dt = w' 2 dt, 

we get that 

f Qxw'dt = (A Ma k+^\A\ 2 k + aha]fd i:j k) [ w' 2 dt + a°-<9,MA I w"'w'dt. 
J~R 2 Jm il 

(11.11) 

Next we will compute L Qiw' dt. We recall that, setting Lq(i/j) — + f'(w)ip, 
the functions ipo and ipi m (1H-7P satisfy 

Lo(i/j )=tw'(t), L (ip 1 )=w". 

Differentiating these equations we get 

Wo) + f"(w)w'fa = (tw'Y, L ($) + - w'" . 

Integrating by parts against w' , using Lq(w') = we obtain 

f"(w)w ,2 i> = - I tw"w' = \J™'\ J f"(w)w' 2 ^ = 
hw'dt = / f'iw^kw' 2 dt = 



Therefore 



a 2 k\A\ 2 / f"{w)^ w' 2 dt - a 2 a%dih d 3 h k / /" '(w)ip lW ' 2 dt = 
Jm Jm 

a 2 k\A\ 2 ]- [ w' 2 - a 2 a\ '.d^djho k [ w"'w'. (11.12) 
2 Jm ' Jm 

Thus, combining relations Ijll.l0|) - (|11.12|) we get 

/ {Qi + --- + Qi)w' dt = (A Ma k + \A\ 2 k + ahalfdijk) f w' 2 dt . (11.13) 
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On the other hand, we observe that 

/ (Q 5 + Q 6 )w'dt = 0(a 2 r- 2 )d ZJ k + O^r- 3 )^ . (11.14) 

J\t+hi\<p a 

Combining relations (|11.13p . (|1 1 . 14)) and (|11.9[) . expansion (jTT75j) follows. Finally, 
for v replaced by ngk w' we have that 

L{kwrjs)rjs kw 1 dt = J rj 2 L(kw)kw dt + J r]s(Aris kw' + 2VrjsV (kw')) kw dt . 

The arguments above apply to obtain the desired expansion for the first integral 
in the right hand side of the above decomposition. The second integral produces 
only smaller order operators in k since Arjs, Vr^ are both of order 0(r~ a ) inside 
their supports. The proof is concluded. □ 

Let us consider now the region 

W := {x e Af 5 I r a (y) < R}, 
where R is a given large number. 

Lemma 11.2. Let k(y) be a smooth function in M a that vanishes when r a (y) = R, 
and set v(y,t) :— r)s(y,t) k(y) w'(t). Then the following estimate holds. 

Q{v,v) = f \Wv\ 2 - f'(u a )v 2 dx = 
Jw 

[ [\V Mc ,k\ 2 ~a 2 \A{ay)\ 2 k 2 } dV a [ w' 2 dt 

+ Ola[ [\Wk\ 2 + a 2 (l + r 4 a )- 1 k 2 ]dV a ) . (11.15) 

V Jr a (y)<R 



Proof. Let us estimate first the quantity 

L(kw') kw' dx . 

Let us express the Euclidean element of volume dx in the coordinates (y,t). Con- 
sider one of the charts ^(y), y S lii of M, I = 1, . . . , N introduced in (|2.6p . which 
induce corresponding charts in M a as in (|2.9p . Then, dropping the index /, we 
compute the element of volume through the change of coordinates 

x - X(y, t) = cT^ay) + (t + h)v{ay) . (11.16) 

Below, as in subsequent computations, the computation of the integral in the entire 
region is performed by localization through smooth partition of unity £i, .. . ,£ m 
subordinated to the covering I = 1, . . . , N of M, namely with the support 

of is contained in Y[(Uk) and = 1- We perform typically a computation of 

an integral of a function g(x) defined in Af$ as 

f g{x)dx = Y J l Zi(ay)g(y,t)dx(y,t) . (11.17) 
J Ms l=1 J Ms 

Let us keep this in mind particularly for estimates obtained by integration by parts 
relative to local variables for y. 
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Let us consider the coordinates X in (|11.16|) . Then we have that 
DX(y,t) = [d x Y + a(t + h)d x v + ad x hv \ d 2 Y + a(t + h) d 2 v + ad 2 h v \ v\ (ay), 
hence 

detDX(y,t) = det [d x Y + a(t + h) d x v \ d 2 Y + a(t + h)h d 2 v \ v] (ay) = 
det[dxY | d 2 Y \ v] + a 2 (t + hf det [d x v \ d 2 v \ v] + 
a(t + h) {det [d x Y | d 2 v | v] + det [d x v \ d 2 Y \ v] } . 
Since mean curvature of M vanishes and the Gauss curvature equals \ A\ 2 , we obtain 

dx = |det DX(y, t) \ dydt = 
(l + a 2 (t + h) 2 \A\ 2 ) |det[9 x y | d 2 Y \ v] | (ay) dydt = 
(l + a 2 (t + h) 2 \A(ay)\ 2 ) dV a (y)dt . 
Using this, we estimate 

I=[ L(kw')kw'dx = 
Jw 

[L(kw')kw'} (l + a 2 (t + h) 2 \A(ay)\ 2 ) dV a (y) dt . 

lr a (y)<RJ\t+h 1 \<p c , 

According to Lemma Til. 11 

7= f(A Ma k + a 2 \A(ay)\ 2 k)k dV a (y) + 



h 

[0( a r- 2 logr a )kdijk + 0( a 2 r^kdjk + 0( a 3 r^k 2 ] dV a (y) 
r a {y)<R 



h 

o ;2 / L(kw') kw' (t + h) 2 \A(ay)\ 2 dV a (y) dt 
w 



Integrating by parts, we see that 

-h= [ [\V Ma k\ 2 -a 2 \A(ay)\ 2 ] dV a (y) . 

The quantity I 2 involves some abuse of notation since it is expressed in local coor- 
dinates for y associated to each chart, and the total should be understood in the 
sense (|11.17[) . Integrating by parts in those coordinates, we get 

h= I [0(ar- 2 logr a )d l kd J k + 0(a 2 r- 3 \ogr a )d l kk + 0(a 3 r-^k 2 } dV a (y) . 

Jr a (y)<R 

Now, 

\a 2 0(r~ 3 \ogr a )d l kk \ < C [a\Vk\ 2 + a 3 (l + r 4 a )" x k 2 } , 
and hence we have 

\h\ < Ca [ [|Vfc| 2 + a 2 (l + ri)- 1 k 2 } dV a (y) (11.18) 

Jr a (y)<R 
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where C is independent of v and all small a. Finally, to deal with the term I3, 
we consider the expression (|11.8|) for L(kw') and integrate by parts once the terms 
involving second derivatives of k. Using that \A\ 2 = 0(r~ 4 ) we then get that 

\h\ < Ca f [ |Vfc| 2 + ^(l + r*)- 1 * 2 ] dV a (y) ■ (H-19) 

Jr a {y)<R 

The same considerations above hold for kw' replaced by rjskw', at this point we 
observe that since nskw' satisfies Dirichlct boundary conditions, we have 

L(kw'r]s) rjgkw' dx = / \V (kw' rjs) \ 2 — a 2 \ri$kw' \ 2 dx = 
iv J w 

{ |VM Q fc| 2 - a 2 \A(ay)\ 2 k 2 } dV a [ w' 2 dt 

r a {y)<R Jm 

+ 0\af [\Wk\ 2 + a 2 (l + r 4 a )- 1 k 2 }dV a ) 

\ Jr a (y)<R J 

and estimate (111 - 15[) has been established. This concludes the proof. □ 

After Lemma Til. 21 the inequality 

m{u a ) > i(M) (11.20) 

for small a follows at once. Indeed, we showed in ^TUIthat the Jacobi operator has 
exactly i(M) linearly independent bounded eigenfunctions Zi associated to negative 
eigenvalues Aj of the weighted problem in entire space M. According to the theory 
developed in we also find that Vz^ = 0(r~ 2 ), hence we may assume 



Q(Zi,Zj) = Xi / qz iZ] dV . (11.21) 



M 



Let us set ki(y) := Zi(ay). According to Lemma fll.21 setting Vi(x) = ki(y) w'(t)r]s 
and changing variables we get 

Q(vi,vj) =a 2 Q(zi, Zj ) fw' 2 + 0(a 3 ) V f |Vz ; | 2 + (1 + r 4 )" 1 ^ 2 dV . (11.22) 
Jr l=i j J M 

From here and relations (jll.21[) . we find that the quadratic form Q is negative on 
the space spanned by the functions V\, . . . , v^ M y The same remains true for the 
functions v t smoothly truncated around r a (y) — _R, for very large R. We have 
proven then inequality (|11.20p . 

In what remains this section we will carry out the proof of the inequality 

m(u a ) < i[M). (11.23) 

Relation (|11.22p suggests that associated to a negative eigenvalue \ of problem 
(|11.3p . there is an eigenvalue of (jll.ip approximated by ~ A^a 2 . We will show next 
that negative eigenvalues of problem (jll.ip cannot exceed a size 0(a 2 ). 

Lemma 11.3. There exists a /1 > independent of R > and all small a such 
that if A is an eigenvalue of problem then 

A > —fi a 2 . 
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Proof. Let us denote 

Qn(^,V) := / |W| 2 - /> Q >/> 2 ■ 
./n 

Then if ip(x) is any function that vanishes for \x'\ > then we have 



Q(i/>,ip) > QjV s n{ rci (y)<R} (4>, i>) + 7 / i> 

3 Ws 



2 



where 7 > is independent of a and R. We want to prove that for some /i>0wc 
have in fl=Afgf] {r a (y) < R} that 

Equivalently, let us consider the eigenvalue problem 



Qfi(V>,V0 > -m 2 I tt — 4 ^ • ( 1L24 ) 



L(V>) + Ap(ax)V> = in Q, (11.25) 
■0 = on r a — R, d n ip = on \t + = p Q . 
Then we need to show that for any eigenfunction ip associated to a negative eigen- 
function, inequality (|11.24|) holds. Here d n denotes normal derivative. Let us 
express this boundary operator in terms of the coordinates {t,y). Let us consider 
the portion of dAfs where 

t + h 1 (ay) = Pa (y). (11.26) 

We recall that for some 7 > 0, p a (y) — p{ay) = 7log(l + r a (y)). Relation (|1 1 .26|) 
is equivalent to 

z - h Q (ay) - Pa (y) = (11.27) 
where z denotes the normal coordinate to M a . Then, for V = V x , we have that a 
normal vector to the boundary at a point satisfying (|11.27p is 

n = Vz - V i; J/!o + p) = v{ay) - aV M (h + p)(ay). 

Now, we have that dtip = V x ip ■ v(ay). Hence, on points (|11.26p . condition d n ip = 
is equivalent to 

d t ^-aWM(ho+ P )-W Ma i> = 0, (11-28) 

and similarly, for 

t + h^ay) =p a (y). (11.29) 

it corresponds to 

d t i)-aV M (ho- p)- V Mc > = ■ (11-30) 
Let us consider a solution tp of problem (|11.25|) . We decompose 

ij) = k(y)w'(t)i] S + V' ± 
where ijs is the cut-off function (|3.12p and 



ip x (y, t) w'(t) dr = for all y E M a n {r a (y) < R}, 

' \r+hi{oty)\<p a {v) 

namely 

Jm w'[tyr]sdt 

Then we have 

Qa(ip, i>) = Qn^,^) + Qu{kw'rjs, kw'rjs) + 2Q n (kw'r ] s, ^ x ). 



56 MANUEL DEL PINO, MICHAL KOWALCZYK, AND JUNCHENG WEI 

Since i/j satisfies the same boundary conditions as ip we have that 
Q a (ilr L ,il> ± ) = - I ( i> x ^ X + f'(u a )ip x2 ) dx . 

Thus, 



Q 



Gntyr 1 .^) - - / / [^ ± A x ^ ± +f(u a )^ ±2 ](l+a 2 (t+h) 2 \A\ 2 )dV a dt. 

Jr a <R J\t+hi\<p a 

Let us fix a smooth function H(t) with H(t) = +1 if t > 1, H{t) = -1 for t < — 1. 
Let us write 

-A^ x - J'K)^ = -9t^ x - f'{w)^ x + ad t [V M (/io + ff(i)p) ■ V^i 1 ] 

Then, integrating by parts in t, using the Neumann boundary condition, we get 
that the integral 

I := 

[dtt^ + /'H^ - ad t (V M (h + H(t)p) ■ V Mq V^)] ^ (l+a 2 (t+h) 2 \A\ 2 ) dt 

\t+ht\<p a 

[d t ^ x - a (V M (ho + H(t) P )-V Ma ip X )] d t iP x (l+a 2 (t+h) 2 \A\ 2 )dt 

\t+h 1 \<p a 

f'{w)iP ±2 (l + a 2 (t + h) 2 \A\ 2 )dt 

\t+hl\<p a 

[d t ^ x - a (V M (h + H{t)p) ■ Vm^ 1 ) ] ^ X 2a 2 (t + h) 2 \A\ 2 dt 

t+fii|<p Q 

[ \d t ilj x \ 2 -f{w)\i) x \ 2 } (l+o(l) ) + aO{r- l )V Ma ^ X d t ^ x +o{l)d t ^ X ^ X dt . 

'\t+hi\<p a 

Now we need to make use of the following fact: there is a 7 > such that if a > is 
a sufficiently large number, then for any smooth function £(t) with J_ a £w' dt = 
we have that 

f ~ f'Me dt > 7 f f + edt. (11.32) 

J — a 

Inequality (| 1 1 .32[) is just a perturbation of the inequality (15. 2| . Wc leave the 
details to the reader. 
Hence 

I > ~ / [\d t ip x \ 2 + \ilj x \ 2 }dt+ f aOir-^Vu^dt^dt. 

(11.33) 

On the other hand, for the remaining part, integrating by parts in the y variable 
the terms that involve two derivatives of ip^- we get that 

II := - [ dt f {A Ma i) X + B^)^{l + a 2 {t + K) 2 \A\ 2 )dV a {y) > 

J\t+hx\<p a Jr a {y)<R 

dt [ |V Mq ^| 2 + o(l)(i> ±2 + \d t i> x \ 2 \V M ^ X \ 2 ) ■ (H.34) 

\t+hi\<p a Jr a (y)<R 
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Using estimates (|11.33|) . (|11.34j) . we finally get 

Qq(VA^) > 3^ / ({dt^l 2 + \V Ma ip ± \ 2 +ip ±2 )dx , (11.35) 
Jn 

for some fi > 0. 

Now, we estimate the crossed term. We have 

-Qa{ilr L ,ku/ris) = [ L(kw'rjs) tp ± (1 + a 2 (t + h) 2 \A\ 2 ) dV a dt . 



Let us consider expression () 1 1 . 8(1 for L(kw'), and let us also consider the fact that 

L(T)skw') = r] S L(kw r ) + 2Vr] S V(kw') + Atj s kw' , 

with the last two terms producing a first order operator in k with exponentially 
small size, at the same time with decay 0(r~ 4 ). Thus all main contributions come 
from the integral 

J} »/.-L /1 i „2/ + i 7 \2 1 a |2 



r} S L(kw') (1 + a 2 (t + h) 2 \A\ 2 ) dV a dt. 

Examining the expression (|11.8p , integrating by parts once in y variable those terms 
involving two derivatives in k, we see that most of the terms obtained produce 
straightforwardly quantities of the type 



o(l) / (\Vk\ 2 + a 2 \A\ 2 k 2 ) dV a + o(l) I (l^l 2 + |VVT). 

J M a 



In fact we have 



I = A Ma kw'risi[> ± dV a dt+ / a 2 a° J d l h d j h kw'" ^dV a dt 
Jn Jn 



h h 
f"(w)4>i kw'i> x dV a dt +6. 



h 

On the other hand, the orthogonality definition of ip 1 - essentially eliminates I\. 
Indeed, 

h = - I A Ma kw' (1-T} S ) il) ± dV a dt= I S7 Ma kw , [{l~r 1 s)\7M a ip ± -ym^)dV a dt = 
Jn Jn 

On the other hand, for a small, fixed number v > we have 

\h\<Ca 2 [ —^r\k\\w"'\\^\dV a dt < Cv-^a 2 [ —-— k 2 dV a +v [ U\ 2 dx . 
Jn l + r l Jmc 1 + r i Jn 

A similar control is valid for I3 since <p-y = 0(a 2 r~ 2 ). We then get 

-1„.2 f 1 7„2 JTA ,. f L/.-1-I2 



I > -Cv-'a* / T k z dV a ~v / |-0 1 • (11.36) 

JM a 1 + r a Jn 

Finally, we recall that from Lemma Til .21 

Qn(kw'r] S ,kw'r)s)= [ [ \V Ma k\ 2 - a 2 \A(ay)\ 2 k 2 } dV a [ w' 2 dt + 9. 

Jr n lv')<R JR 



(11.37) 
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From estimates (|11.35|) . (|11.36[) . (|11.37|> . we obtain that if v is chosen sufficiently 
small, then 

Qn(^,^)>-Ca 2 f -^—k 2 dV a >-/i« 2 f -J—\^dx, 

JM a 1 + r a Jn 1 + r a 

for some /x > and inequality (|11.24jl follows. □ 

In the next result, we show that an eigenfunction with negative eigenvalue of 
problem (jll.ip or l|11.2[) decays exponentially, away from the interface of u a . 

Lemma 11.4. let <p be a solution of either (| 1 1. 1[ or (1 1 1 . 2\ with A < 0. Then <j> 
satisfies in the subregion of Af a where it is defined that 

\<Kv,t)\ < C e- aW (11-38) 

where a > can be taken arbitrarily close to min{<7+, cr_}. The number C depends 
on a but it is independent of small a and large R. We have, moreover, that for 
\ax'\ > Rq, 

rn 

\(f)( x )\<C^e- alx3 - a ~ 1{Fk{ax ' )+f3 J alos ^ x ' l)l . (11.39) 
i=i 

where Rq is independent of a. Finally, we have that 

\Mx)\<Ce- a i for dist (x, M a ) > — . (11.40) 

a 

Proof. Let solve problem (jll.ip for a large R. Let us consider the region between 
two consecutive ends Mj >a and Mj+i, Q . For definitcncss, we assume that this region 
lies inside S+ so that f'(u a ) approaches cr 2 . inside it. So, let us consider the region 
S of points x — (x',X3) such that r a {x) > Rq for a sufficiently large but fixed 
Ro > and 

{a.j 7 a(3j) \oga\x'\ 7 bj 7 cry < ax 3 < (a j+ i 7 a(3. ]+ i ) loga|a;'| + b J+ i — cry. 

In terms of the coordinate t near -Mj jQ , saying that 

ax3 ~ (aj 7 a(3j) log a\x \ +bj 7 cry 

is up to lower order terms, the same as saying t ~ 7, similarly near Mj+i jtx . Thus 
given any small number r > we can choose 7 sufficiently large but fixed, inde- 
pendently of all Rq sufficiently large and any small a, such that 

f(u a ) < -{o-+-Tf in 5. 

Let us consider, for x £ S and a = 07 — 2r the function 

V\{x) :— e - a [ x 3- a ~ 1 ( a i+ a 0j) lo & a \ x '\+ b 3] 7 e -cr(a- 1 [a j + 1 +a/3 j + 1 )lo g a\x'\+b j + 1 )-x 3 ] 

Then v has the form 

wi = A x e- ax *r A * 7 B x e aX3 r- B \ r = \x'\, 

so that 

Avi = Ajr- 2 r A2 A ie - aX3 7 Bfr" 2 B ir - B2 e aX3 7 a 2 v x < 
[a 2 A 2 2 R„ 2 + a 2 B 2 Ro 2 + a 2 ] v 1 . 

Here 

A 2 = aa~ 1 (aj 7 a/3j), B 2 = aa~ (dj+i 7 a(3j + i). 
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Hence, enlarging Rq if necessary, we achieve 

Avi + f'(u a )vi < in S. 
Therefore v so chosen is a positive supersolution of 

Av + f'(u a )v + Xp(ax)v < in S. (H-41) 
Observe that the definition of v also achieves that 

inf > 7 > 

dS\{r a =R } 

where 7 is independent of a. Now, let us observe that the function v 2 = e~"" x >~~> 
also satisfies, for small a, inequality (|11.41j) . As a conclusion, for <f>, solution of 
(jll.ip . we have that 

10(a) I < C\\(f>\\ 00 [v 1 (x)+v 2 (x)} forall x e S,r a (x) < R. (11.42) 

Using the form of this barrier, we then obtain the validity of estimate (jll.39[) , in 
particular that of (|11.38[) , in the subregion of Ms in the positive t direction of Mj a 
and Mj+i tOI when r a (y) > Rq- The remaining subregions of Ms H {r a (y) > Ro} are 
dealt with in a similar manner. Finally, to prove the desired estimate for r a {y) < ro 
we consider the region where \t\ < — assuming that the local coordinates are well 
defined there. In this case we use, for instance in the region 

28 

v < t < — 

a 

for v > large and fixed, a barrier of the form 

v(y,t) = e-° t +e- a ^- t K 
It is easily seen that for small a this function indeed satisfies 

A x v - f(u a )v < 
where a can be taken arbitrarily close to cr+. We conclude that 

\4>{y,t)\ < CUW^e-^ for^<t< - . 

a 

Thus estimate (jll.38|) holds true. Inequality () 1 1 .40(1 follows from maximum prin- 
ciple. 

Finally, for a solution of problem (|11.2p the same procedure works, with only mi- 
nor difference introduced. Estimate (|11.42[) can be obtained after adding a growing 
barrier. Indeed, we obtain 

\<f>(x)\ < C\\<f>\\oo [vi(x) +v 2 (x) + sv 3 (x)] forall xeS 

with V3(x) = se a \ x I, and then we let e — ► 0. We should also use ee' 7X3 to deal 
with the region above the last end M m and similarly below Mi. We then use the 
controls far away to deal with the comparisons at the second step. The proof is 
concluded. □ 
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II. 1. The proof of inequality (|10.23|h Let us assume by contradiction that 
there is a sequence a = an — ► along which 

m(u a ) > i(M) =: N . 

This implies that for some sequence R n — ► +00 we have that, for all R > R ni 
Problem (|11.1|) has at least N + 1 linearly independent eigenfunctions 

<t>l,a,R, ■ ■ ■ , 4>N+l,a,R 

associated to negative eigenvalues 

M,a,R < M,a,R < • • • < Aat+1,q,_R < . 

We may assume that ||0i jQ ,.R ||oo = 1 and that 

p(ax)4> iiatR (f> j:atR dx = for all i, j = 1, . . . , N + 1, i £ j. 

Let us observe that then the estimates in Lemma \l 1 . 41 imply that the contribution 
to the above integrals of the region outside A/5 is small. We have at most 

/ p{ay)(j)i >a ^<P],a,Bdx = 0{a 3 ) for all i,j = l,...,N + l, i^j. (11.43) 

From the variational characterization of the eigenvalues, we may also assume 
that \i, a ,R defines a decreasing function of R. On the other hand, from Lemma 

III. 3l we know that Xi ia .R = 0(a 2 ), uniformly in i?, so that we write for convenience 

K,aM = Hi,a,R.O?, (M,a,R < 0. 

We may assume Hi, a ,R. —> < as R — > +00. We will prove that <fii,a,R. 
converges, up to subsequences, uniformly over compacts to a nonzero bounded 
limit 4>i,u which is an eigenfunction with eigenvalue /ii iQ a 2 of Problem (|11.2|1 . We 
will then take limits when a — > and find a contradiction with the fact that J has 
at i(M) negative eigenvalues. 

We fix an index i and consider the corresponding pair <f>i^ a ,R, (J-i,a,R, to which 
temporarily we drop the subscripts i, a, R. 

Note that by maximum principle, \<p\ can have values that stay away from zero 
only inside A/5 Besides, from Lemma Til .41 = 0(e~ <T ' t ') in A/5. We observe then 
that since A remains bounded, local elliptic estimates imply the stronger assertion 

\D 2 <j)\ + \Dcj)\ + \<j)\ < Ce- CT|t| mtf s . (11.44) 

In particular, considering its dependence in R, <fr approaches up to subsequences, 
locally uniformly in R 3 a limit. We will prove by suitable estimates that that 
limit is nonzero. Moreover, we will show that <j) ~ z(ay)w' (t) in A/5 where z is an 
eigenfunction with negative eigenvalue w \x of the Jacobi operator J . 

First, let us localize <p inside A/5. Let us consider the cut-off function r/s in (|3.12p . 
and the function 

4> = m<P- 

Then 4> satisfies 



L{4>) + fia 2 q(ax)(f> = E, 



= -2Vr, s V<t> - A?750 



(11.45) 
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with L(4>) = Acf> + f'(u a )(j). Then from ()L1.44|) we have that for some a > 0, 

\E a \ <Ca 3 e- CT l*l(l+^)- 1 . 
Inside Afs we write in (y, t) coordinates equation (|1 1 .45|) as 

L*(4>) + B(4>) + X P (ay)4> = E a (11.46) 

where 

L*(4>) = d tt 4> + A M J + f(w(t))4> . 

Extending (f> and E a as zero, we can regard equation (|11.46|) as the solution of a 
problem in entire M a x R for an operator L that interpolates L inside Af$ with L* 
outside. More precisely <fi satisfies 

L(4>) := L*{4>) + B{4>) + \p(ay)4> = E a in x R, (11.47) 

where for a function tp(y,t) we denote 

B(ih) ■= ( xB ^ ii \t + hx(ay)\ < p a (y) + 3 (1148) 
' | otherwise 

and 

x(y,t) = {i(y+(t + h)v a (y)) 
with Ci the cut off function defined by (|4.8p for n = 1. In particular, L = L in A/5. 
Now, we decompose 

t) = (p(y,t) + k{y)T] S w'(t) (11.49) 

where 

Ky) = - W \t)kp^± 

J R r) S w> dr 

so that 

/ <p(y,t)w'{t)dt = for all y£M*. 
Jr 

From (|11.44p , k is a bounded function, of class C 2 defined on with first and 
second derivatives uniformly bounded independently of large R. A posteriori we 
expect that k has also bounded smoothness as a function of ay, which means in 
particular that Dk = 0(a). We will see that this is indeed the case. 

The function (p satisfies the equation 

L(<p) + [ia 2 p(ay) ip = — L(kw') + E a — \xa 2 pkw' inAff x R . (11.50) 

We observe that the expansion (|11.8|) holds true globally in x R for L(kw') re- 
placing L(kw'). We also have the validity of expansion (| 11 . 5|) for the corresponding 
projection, namely 



L(kw')w'dt = (A Mc ,k + a 2 \A\ 2 k) [ w' 

Jr 



2 dt 



+ 0(ar- 2 ) dijk + 0(a 2 r- 3 ) d t k + 0{a 3 r - 4 ) k . (11.51) 
Thus, integrating equation (| 1 1 . 50|) against w' we find that k satisfies 

A Ma k + a 2 \A\ 2 k + p,a 2 p(ay)k + 

0(ar- 2 )d tJ k + 0(a 2 r- 3 )d,k + 0{a 3 r~ 4 )k = 

0(a 3 r- 4 ) - —Lj f Bfa) w' dt, y E M* . (11.52) 

Jr w 



(i2 
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Let us consider the function z(y) denned in M by the relation k{y) = z(ay). Then 
(|11.52|) translates in terms of z as 

A M z + \A(y)\ 2 z + f*q(y)z = 

a [0(r- 2 ) d ijZ + 0(r- 3 ) d iZ + 0(r- 4 ) z + 0(r- 4 )] +B y e M R . (11.53) 
where 



B(y) :-- 



1 



Ir w ' 



B(tp)(a~ 1 y,t)w' dt, y E M 1 



(11.54) 



In other words we have that k(y) = z(ay), where z solves "a perturbation" of the 
eigenvalue equation for the Jacobi operator that we treated in 101 We need to 
make this assertion precise, the basic element being to prove that the operator B[z] 
is "small" . For this we will derive estimates for (p from equation (|11.50p . 

We shall refer to the decomposition Qi + ■ ■ ■ + Qe in (|11.8|) to identify different 
terms in L(kw'). Let us consider the decomposition 

where (pi solves the linear problem for the operator L* and the part of L(kw') that 
"does not contribute to projections" , namely 

Qz + Qi = —w" aa^Adjhdik + dihdjk) + a 2 kAMhi + a 2 ha\f [djhdik + dihdjk) 



atw 



<Zjj dijk + ab]'°dik 



(11.55) 



More precisely tpi solves the equation 



L*(<pi) + a 2 fj,pip! = Q 3 + Q4 hxM^xR. (11.56) 
This problem can indeed be solved: according to the linear theory developed, there 
exists a unique solution to the problem 

L*(ipi) + fia 2 p ipi = Q 3 + Qi + c(y)w' (t) in M% x M, 

such that 

/ ipi w' dt = for all y € M% 
Jr 

and 

PV|| Pj l,<7 + ||^l||oo,l,a+ ||Vl||oo,l,«r < ||Qs + Q4pA,* < Ca ■ ( 1L57 ) 

But since 

/ (Qs + Qi) w'dt = for all y e 

JR 

it follows that actually c(y) = 0, namely ipi solves equation (lll.56p . 

We claim that <p% has actually a smaller size than tp±. Indeed tf2 solves the 
equation 

L*(p 2 )+B(ip 2 )+fia 2 p<p2 = E a -B(ip 1 )-(Q 1 +Q 2 +Q 5 +Qe)~^a 2 qkw' in M*xR. 

(11.58) 

Now, we have that 

Qi + Qi + Qa + Qe = 



A.M a k + aha^ d^k 



w' + a 2 [-\A\ 2 ktw" + a^dihodjhokw'" + 
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2 f"(w) 0i kw' + (t + hfaj^jkw' + 2{t + h)a% dihdjk) w" ] + 

0(e^r^ 2 )d tJ k + Oie-^r-^dM = (11.59) 



a 

0(a 2 r~ 2 \og 2 r a e-^) + p(y)w'(t) , 
for a certain function p{y). On the other hand, let us recall that 
B = (f'(u a ) - f'(w)) - a 2 [(t + ti^Af + A M /ii]& -aa%{djhd it + dihd jt ) + 

a(t + h) [aljdij - 2a ajjdihdjt + a(b]di - ab]dihd t ) ) ] + 

a 3 {t + h) 2 b\d t + a 2 [al + a(t + h)a 1 lJ )}d t hd J hd t t (11.60) 

Thus the order of B(y>i) carries both an extra a and an extra r" 1 over those of ipi, 
in the sense that 

||B(¥>i)|| p ,2,<r <Ca 2 . (11.61) 
From relations (|1 1 .59|) and (|11.61[) we find that ip 2 satisfies an equation of the form 

L*(ip 2 ) +B(tp 2 ) + pa 2 qip 2 =g + c(y)w' in M* x R (11.62) 

where for arbitrarily small a' > we have 

||S||P,2-<7',<T < CO 2 . 

Since if2 satisfies L ifi2 w' dt = 0, the linear theory for the operator L» yields then 
that 

\\D 2 <p2\\ p ,2-*>,<, + || D<p2\\oc,2-*>,a + || H|oc,2- CT ', CT < Co? , (11.63) 

which compared with (|11.57j) gives us the claimed extra smallness: 

||B(pa)||p,3-«r',«r < Co?. (11.64) 

Let us decompose in (|11.54l) 

B = B x + B 2 

where 

Bi := -r^a- 2 [ B(ipi){a- 1 y,t)w' dt, 1 = 1,2. (11.65) 
Jr w J* 
From Lemma 18.11 we get that 

||Sl|| p ,2-2- CT < < ^«- 2 ||B(^i)||p j2 ,a < C (11.66) 

and 

l|Ba|| p ,3-| -2.' < Ca- 2 ||B(^ 2 )|| p , 3 _ CT ,, ff < Ca (11.67) 

Now, we apply the estimate in part (b) of Lemma 110.21 to equation (|11.53[) and 
then get for z(y) = fc(^) the estimate 

ll^llp,2-l-2 CT ' + II (1 + \x\f- 2 "' DzU < C [ ||/|| p , 2 -2-2 CT < + Iklloo ] (11-68) 

where 

/ = a [0(r- 2 ) d ijZ + 0(r' 3 ) 8 t z + 0{ r - 4 ) z + 0(r~ 4 )] + B. 
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Then from estimate (| 1 1 .68[) it follows that for small a, 

\\D 2 z\\ p , 2 _2_ 2a , + || (1 + \x\) 1 ^' Dz \\oo < Ca. (11.69) 

Using this new information, let us go back to equation (|11.56[) and to the expression 
(|11.55[) for Q3 + Q4. The terms contributing the largest sizes in this function can 
be bounded by 



Now, we compute 



\Dk\ \D 2 k\ 



l + r a 1 + r 



(1 + ^) 2 ) P / ^$ydV a < 

JB(yS) ( l + r ar 

Ca 2p - 2 f \D 2 z\ p dV < Ca 2p - 2 \\D 2 z\\ p2 _2_ 2a , < Ca 2p ~ 2 , 

JB(y, a ) 

and 

(l + r Q (y) 2 - 2 °'y [ 7^ Vp dV a < 
Jb( v ,i) (1 +r a ) p 

C|||Cfc|(l + r Q ) 1 - 2CT 'f 00 = Co? || [-D*| (1 + r) 1 - 2CT 'f oc < Co? . 
As a conclusion, from expression (jll.55|) we obtain that 

IIQ3 + Q4IU2-2CT' ,<r < Ca 2 , 

and therefore a substantial reduction of the size of ipi, compared with (|11.57[) . we 
have 

PVllp^aV + ||^i||oc,2- 2ct > < Ca 2 , (11.70) 
hence, using again Lemma 18. II we get 

IIBillp.a-a-a^ < Car 2 \\B{<pi)\\ p ,a-2*>,* < Ca (11.71) 

which matches the size we initially found for B 2 in (|11.67[) . 

We recall that <f> — <fii ya ,R has a uniform C 1 bound (|11.44p . Thus, passing to a 
subsequence if necessary, we may assume that 

4>i,cc,R -> <t>i, a as R -> +00, 

locally uniformly, where <^j jQ , is bounded and solves 

&<t>i,a + f'{n a )<j>i, a + M»,a a 2 p(ax) <£ i)Ct = in R 3 . (11.72) 

Let us return to equation (|11.53[) including the omitted subscripts. Thus k = k^ at R 
satisfies the local uniform convergence in M a , 

h.a,R(y) = C 4>i,a,Rw' dt — > C 4>i, a w' dt =: fei, a (y) ■ 

Jlt+/lll<p a Jlt+/lll<p Q 
We have that z = satisfies 

AA/Zi. a ,ji + |A(y)| 2 ^ >a)iJ + Hi, a ,Rq(y) Zi t0lt R = 

a [0(r~ 2 ) r;,„ •„,;,. + 0(r- 3 ) + 0(r~ 4 ) ^ jQ>iJ + 0(r- 4 )] + 

Bi, a<R , yeM R , 

where 

||B ijCe , H || Pi 3_ f _ 3CT , < Ca (11.73) 
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with arbitrarily small a' > and C independent of R. We apply now the estimates 
in Lemma 110.21 for some 1 < p < 2 and find that for C independent of R we have 

\\Zi, a ,R\\L^(M R ) < C [ ||«f,a,ii||ioo( r<fla ) +0(a)] 

or equivalcntly 

||**,a,B|U-.(jif*) <C[\\k i!a>R \\ L «, {ra<Ro) +0(a)} . (11.74) 
Since from l|11.49|) we have that 

<t>i,R,a(v,t) = ipi,R, a (v,t) + h iR . a (y) w' (t) in Afg, r a (y) < R, (11.75) 
where we have uniformly in R 

\<Pi,RAv>t)\=0(ae-' r Mr- 2 ), 
while <j>% R t a = 0(e~~) outside Afg, and 

||0i,fl,a||oo = 1) 

then 

||fci,a,fl||z,»(Af«) > 7 > 
uniformly in R. Thus from (|11.74|) . the limit fcj iQ as R — ► +oo cannot be zero. We 
have thus found that </) !iQ is non-zero. Moreover, we observe the following: Since 
the functions 

Zi := diU a , i = 1,2,3, Z 4 := -x 2 d\u a + xid 2 u a 
are bounded solutions of (| 11 . 2[) for A = 0, we necessarily have that 

/ p(ax) Zj4>i, a dx = 0, j = 1,2,3,4. 

JR 3 

Let 

4 

= 22 dqZi, i = 1, . . . , J. 
i=i 

Then we also have 

/ p{ax) Zi4>i, a dx — 0, i=l,...,J. (11.77) 

Jm 3 

Now we want to let a — ► 0. satisfies 

&MZi,a + \A(y)\ 2 Z iia + fii,aP(y) Zi,a = 

a [0(r" 2 ) %^ i>a + 0(r- 3 ) d t z ua + 0(r" 4 ) z ijQ + 0(r~ 4 )] + B i>a , y £ M , 
with 

l|Bi,a|| p , 3 - = - 3ff < < Ca. (11.78) 

Moreover, 

IN,a||i>»(Af) < C [ ll z i,a,fllU=°(r<_R, ) + O(a) ] . 

Since we also have that 

H-D 2z i,a||LJ>(A/) < C i 

Sobolev's embedding implies that passing to a subsequence in a, Zi tCt converges as 
a — > 0, uniformly over compact subsets of M to a non-zero bounded solution Zi of 
the equation 

A M Zi + \A{y)\ 2 z t + mq(y)zi = in M, 

with Hi < 0. 



(11.76) 
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Now, we have that 



where 

We recall that 

Since on Ms , 
we get then that 

or 



\Vi(y,t)\ <Caer*W. 
q(ay)<pi a (f>j a dx = O(a) for all i ^= j 
dx = (l + a 2 \A\ 2 (t + h))dV a dt, 



q(ay)z it a(ay) z ha (ay) dV a = 0(a) 



q{y)zi, a {y) Zj, a (y) dV = 0(a ) for all i ^ j. 

I M 

We conclude, passing to the limit, that the z,'s i = 1, . . . , N + 1 satisfy 

q ZiZj dV — for all i =/= j. 

M 

Since, as we have seen in qiQ[ this problem has exactly N = i(M) negative eigen- 
values, it follows that fJ.N+1 = 0, so that that z^v+i is a bounded Jacobi field. 

But we recall that, also 

Z t = Zi (ay)w'(t) +0{ae-^ tl ) for all i = 1, . . . , J, 

hence the orthogonality relations (|11.77p pass to the limit to yield 

q z,i ■ zjv+i dV — 0, i — 1, . . . , J . 

M 

where Zi's are the J linearly independent Jacobi fields. We have thus reached a 
contradiction with the non-degeneracy assumption for M and the proof of m(u a ) — 
i(M) is concluded. 

Finally, the proof of the non-degeneracy of u a for all small a goes along the same 
lines. Indeed, the above arguments are also valid for a bounded eigenfunction in 
entire space, in particular for /i — 0. If we assume that a bounded solution Z5 of 
equation ()11.2j) is present, linearly independent from Zx, . . . , Z4, then we assume 
that 

/ p(ax)Z 5 Zidx = i = l,...,J. (11.79) 

Thus, if in the same way as before, we have that in J\f$, 

Z 5 = z 5 (ay)w'(t) + tp 

with tp orthogonal to w' (t) for all y and ip small with size a and uniform exponential 
decay in t. The function z a solves an equation of the form (|11.78p . now for /1 = 0. 
In the same way as we did before, it converges uniformly on compacts to a non-zero 
limit which is a bounded Jacobi field. But the orthogonality (|11.79|1 passes to the 
limit, thus implying the existence of at least J + 1 linearly independent Jacobi field. 
We have reached a contradiction that finishes the proof of Theorem [21 □ 
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12. Further comments and open questions 

12.1. Symmetries. As it is natural, the invariances of the surface are at the same 
time inherited from the construction. If M is a catenoid, revolved around the X3 
axis, the solution in Theorem [T] is radial in the first two variables, 

u a (x) = u a ( \x'\, x 3 ) . 

This is a consequence of the construction. The invariance of the Laplacian under 
rotations and the autonomous character of the nonlinearity imply that the entire 
proof can be carried out in spaces of functions with this radial symmetry. More 
generally, if M is invariant a group of linear isometries, so will be the solution found, 
at least in the case that f(u) is odd. This assumption allows for odd reflections. 
The Costa-Hoffmann-Meeks surface is invariant under a discrete group constituted 
of combination of dihedral symmetries and reflections to which this remark apply. 

12.2. Towards a classification of finite Morse index solutions. 

Understanding bounded, entire solutions of nonlinear elliptic equations in l w is a 
problem that has always been at the center of PDE research. This is the context of 
various classical results in PDE literature like the Gidas-Ni-Nirenberg theorems on 
radial symmetry of one-signed solutions, Liouville type theorems, or the achieve- 
ments around De Giorgi conjecture. In those results, the geometry of level sets of 
the solutions turns out to be a posteriori very simple (planes or spheres). More 
challenging seems the problem of classifying solutions with finite Morse index, in a 
model as simple as the Allen-Cahn equation. While the solutions predicted by The- 
orem [T] are generated in an asymptotic setting, it seems plausible that they contain 
germs of generality, in view of parallel facts in the theory of minimal surfaces. In 
particular we believe that the following two statements hold true for a a bounded 
solution u to equation (jl.lj) in R 3 . 

(1) If u has finite Morse index and V«(s) 7^ outside a bounded set, then each 
level set of u must have outside a large ball a finite number of components, each 
of them asymptotic to either a plane or to a catenoid. After a rotation of the 
coordinate system, all these components are graphs of functions of the same two 
variables. 

(2) If u has Morse index equal to one. Then u must be axially symmetric, namely 
after a rotation and a translation, u is radially symmetric in two of its variables. 
Its level sets have two ends, both of them catenoidal. 

It is worth mentioning that a balancing formula for the "ends" of level sets to the 
Allen-Cahn equation is available in R 2 , see [15]. An extension of such a formula to 
R 3 should involve the configuration (1) as its basis. The condition of finite Morse 
index can probably be replaced by the energy growth (|1.9[) . 

On the other hand, (1) should not hold if the condition V11 7^ outside a large 
ball is violated. For instance, let us consider the octant {x\,X2,x^ > 0} and the 
odd nonlinearity f(u) = (1 — u 2 )u. Problem (jl.lj) in the octant with zero boundary 
data can be solved by a super-subsolution scheme (similar to that in [5]) yielding a 
positive solution. Extending by successive odd reflections to the remaining octants, 
one generates an entire solution (likely to have finite Morse index), whose zero 
level set does not have the characteristics above: the condition Vit 7^ far away 
corresponds to embeddedness of the ends. 
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Various rather general conditions on a minimal surface imply that it is a catenoid. 
For example, R. Schoen |35j proved that a complete embedded minimal surface in 
R with two ends must be catenoid (and hence it has index one) . One may wonder 
if a bounded solution to (jl.ip whose zero level set has only two ends is radially 
symmetric in two variables. On the other hand a one-end minimal surface is forced 
to be a plane [20] . We may wonder whether or not the zero level set lies on a half 
space implies that the solution depends on only one variable. 

These questions seem rather natural generalizations of that by De Giorgi, now on 
the classification finite Morse index entire solutions of Ijl.ljl . The case in which the 
minimal surfaces have finite topology but infinite total curvature, like the helicoid, 
are natural objects to be considered. While results parallel to that in Theorem [T] 
may be expected possible, they may have rather different nature. The condition of 
diverging ends in f3 is not just technical. If it fails a solution may still be associated 
to the manifold but interactions between neighboring interfaces, which are inherent 
to the Allen-Cahn equation but not to the minimal surface problem, will come into 
play. The case of infinite topology may also give rise to very complicated patterns, 
we refer to Pacard and Hauswirth [17j and references therein for recent result on 
construction of minimal surfaces in this scenario. 

13. Appendix 

In this appendix we carry out the computations that lead to Lemma 12.11 

13.1. Coordinates near M and the Euclidean Laplacian. Let us consider the 
smooth map 

(y, z) e M X R i — > x = X(y, z)=y + zv(y) e R 3 . (13.1) 
Let O be a set as in the statement of Lemma |2~T1 and consider the subset of M x R 
defined as 

6 = {{ay,a{t + h{y)))e M xR/ {t,y)eO}. 
Then X\q is one to one, and 

6 c {(y, z) e M x R / \z\ < 6 log(l + r(y))}. 

Since along ends diV = 0(r~ 2 ) so that zdiV is uniformly small in O, it follows that 
X is actually a diffeomorphism onto is image, TV = X(0) = aj\f. 

The Euclidean Laplacian can be computed in such a region by the well-known 
formula in terms of the coordinates [y, z) e O as 

A, = d zz + A Mz - fljw,0». x = *(tf> 2 )> (V> z )^° ( 13 - 2 ) 
where M z is the manifold 

M z = {y + zv{y) /ye M}. 

Local coordinates y = ifc(y), y € R 2 as in (12. ip induce natural local coordinates in 
M z . The metric gij(z) in M z can then be computed as 

gij (z) = (diY, djY) + z«a i F, 8 jV ) + (djY, dp)) + z 1 {diu, d jV ) (13.3) 

or 

gij (z)= gij +zO(r- 2 ) +z 2 0(r- i ). 
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where these relations can be differentiated. Thus we find from the expression of 
Am, in local coordinates that 

A Mz = A M + za^iy, z)d ij + zb\{y, z)d u y = Y(y) (ISA) 

where a\j,b\ are smooth functions of their arguments. Let us examine this expan- 
sion closer around the ends of Mj, where y = Y^(y) is chosen as in (I2.1|) . In this 
case, from (| 13 . 3[) and (|2.2p we find 

gV(z) = g v +zO(r- 2 ) + z 2 0(r 4 ) + ... 

Then we find that for large r, 

A Ma = A M + zO( r - 2 )d lJ + zO(r- 3 )d t . (13.5) 

Let us consider the remaining term in the expression for the Laplacian, the mean 
curvature Hm^ ■ We have the validity of the formula 

H M ;=± T ^- z =±k l + k 2 z + klz 2 + ... 

i=l 1 i=l 

where fc^, i = 1,2 are the principal curvatures. Since M is a minimal surface, we 
have that k\ + \i2 = 0. Thus 

|A| 2 = k\ + k\ = -2kik 2 = -2K 

where \A\ is the Euclidean norm of the second fundamental form, and K the Gauss 
curvature. As r — > +oo we have seen that fej = 0(r~ 2 ) and hence \A\ 2 = 0(r~ A ). 
More precisely, we find for large r, 

H Mz = \A\ 2 z + z 2 0{r- 6 ). 

Thus we have found the following expansion for the Euclidean Laplacian, 

A x = d zz + A M - z\A\ 2 d z + B (13.6) 

where expressed in local coordinates in M the operator B has the form 

B = za^y, z)dij + z b\{y, z)d l + z 2 b\{y, z)d z (13.7) 

with ajj, bj, b\ smooth functions. Besides, we find that 

a\ J {y : z) = 0{r- 2 ), b](y, z) = 0{ r -% bl(y,z) = 0( r - 6 ), (13.8) 

uniformly in z for (y, z) G O. Moreover, the way these coefficients are produced 
from the metric yields for instance that 

a lj(yy z ) = a lj(y, o) + za [ 2 j (y> z )> a h{y> z ) = o{r~ 3 ), 

bl(y,z) = bl(y,0) + zb^(y,zl bf\y,z) = 0^) . 
We summarize the discussion above. Let us consider the parameterization in 
(EDO]! of the region N. 

Lemma 13.1. The Euclidean Laplacian can be expanded in Af as 

A x = d zz + — Hm^9 z = 
d zz + A M - z \A\ 2 d z + z [4(y, z)d ij + &J(y, z)d l ] + z 2 b\{y, z)d z , 
A M = af-jdij + b° i d ll x = X(y, z), (y, z) eO, 
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where Oy , bj are smooth, bounded functions, with the index k omitted. In addition, 
for k — 1, . . . , m, 

a\ 3 = S ZJ d Q i + 0{r-% b\ = 0(r- 3 ), b\ = 0{r- & ) , 

as r — |y | — > oo, uniformly in z variable. 

13.2. Laplacian in expanded variables. Now we consider the expanded minimal 
surface M a = ar x M for a small number a. We have that Af = a~ 1 Af. We describe 
Af via the coordinates 

x=X(y,z):=y + zu a {y), {y,z)ea- 1 6. (13.9) 

Let us observe that 

X(y, z) — a~ 1 X(ay, az) 

where x = X(y,z) = y + zv(y), where the coordinates in Af$ previously dealt with. 
We want to compute the Euclidean Laplacian in these coordinates associated to 
M a . Observe that 

A x [u(x)} \ x = X (y,z) = a 2 Ai[u(a _1 i)] \^ x{ay>az) 

and that the term in the right hand side is the one we have already computed. In 
fact setting v(y, z) := u{y + zv a (y)), we get 

Axu\ x=X (y,z) = a 2 {Ay,Mz+dzz-H Ms d s ) [v(a' x y, a' 1 z)] | (fi)2 ) =(aj , jCw) • (13.10) 

We can then use the discussion summarized in Lemma ll3.1l to obtain a represen- 
tation of A x in Af via the coordinates X(y,t) in (|13.9p . Let us consider the local 
coordinates Yk a of M a in (|2.9|) . . 

Lemma 13.2. In Af we have 

A x = d zz + &M a ,x ~ HM a:Z d z = 

d zz +AM a —a 2 z | A{ay)\ 2 d z + az [a}j (ay, az)dij +ab] (ay , az)di] +a 3 z 2 b\(ay, az)d z , 

A M a = a° lJ (ay)d lj + b\(ay)d ll (y, z) g a^O, y = Y ka (y) 
where a-, b 1 ^ are smooth, bounded functions. In addition, for k — 1, . . . ,m, 

a\. = 8 l3 5 a + 0(r- 2 ), b\ = 0{r~% b\ = 0(r" 6 ) , 
QS r a{y) = \ocy\ — > oo, uniformly in z variable. 
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13.3. The proof of Lemma 12.11 Let us consider a function u defined in Af, 
expressed in coordinates x — X(y,z), and consider the expression of u in the 
coordinates x = Xh{y,t), namely the function v(j,t) defined by the relation in 
local coordinates y = Yfc(y), 

v(y,z - h(ay)) =u[j,z), 

(by slight abuse of notation we are denoting just by h the function ho Y" fe ). Then 
we compute 

d L u = diV — adtvdih, d z u = dtV, 
dijU = dijV — aduvdjh — adjtvdih + a 2 duvdihdjh — a 2 dtvdijh . 
Observe that, in the notation for coefficients in Lemma ri3.2[ 

a°jdijh + b®d t h = A M h, + ab^div = A Ma v ■ 

We find then 

A, = da + A Ma - a 2 [(t + h)\A\ 2 + A M h]d t - 2a a% d 3 hd lt + 
a(t + h) [a]jdij — 2aajj dihdjt + a{b\di — ab\dihdt) ] + 

a 3 (t + h) 2 b\d t + a 2 [a°- + a(t + ]d t hd 3 hd tt (13.11) 

where all the coefficients are understood to be evaluated at ay or (ay, a(t + h(aj)) . 
The desired properties of the coefficients have already been established. The proof 
of Lemma 12.11 is concluded. □ 
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